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ABSTRACT 

We present a re-analysis of the XMM-Newton long-look of the X-ray bright Seyfert galaxy IC 4329a. 
The Fe K bandpass is dominated by two peaks, consistent with emission from neutral or near-neutral 
Fe Ka and K/3. A relativistic diskline model whereby both peaks are the result of one doubly-peaked 
diskline profile is found to be a poor description of the data. Models using two relativistic disklines 
are found to describe the emission profile well. A low-inclination, moderately-relativistic dual-diskline 
model is possible if the contribution from narrow components, due to distant material, is small or 
absent. A high-inclination, moderately relativistic profile for each peak is possible if there are roughly 
equal contributions from both the broad and narrow components. Upper limits on Fe XXV and Fe XXVI 
emission and absorption at the systemic velocity of IC 4329a are obtained. We also present the results 
of RXTE monitoring of this source obtained so far; the combined XMM-Newton and RXTE data sets 
allow us to explore the time-resolved spectral behavior of this source on time scales ranging from hours 
to 2 years. We find no strong evidence for variability of the Fe Ka emission line on any time scale 
probed, likely due to the minimal level of continuum variability. We detect a narrow absorption line, 
at a energy of 7.68 keV in the rest frame of the source; its significance has been confirmed using Monte 
Carlo simulations. This feature is most likely due to absorption from Fe XXVI blueshifted to ~0,lc 
relative to the systemic velocity, making IC 4329a the lowest-redshift AGN known with a high-velocity, 
highly-ionized outflow component. As is often the case with similar outflows seen in high-luminosity 
quasars, the estimated mass outflow rate is larger than the inflow accretion rate, signaling that the 
outflow represents a substantial portion of the total energy budget of the AGN. The outflow could arise 
from a radiatively-driven disk wind, or it may be in the form of a discrete, transient blob of ejected 
material. 

Subject headings: galaxies: active — galaxies: Seyfert — galaxies: individual (IC 4329a) 


1. INTRODUCTION 

The hard X-ray emission of Seyfert 1 AGN is dominated 
by rapidly- variable emission thought to originate from in- 
verse Comptonization of soft seed photons by a hot corona 
near the central black hole (e.g., Shapiro, Eardley & Light- 
man 1976; Sunyaev & Titarchuk 1980; Haardt et al. 1994). 
Furthermore, the disk, or some other cold, optically thick 
material, reprocesses the hard X-rays, as evidenced by the 
so-called ’Compton reflection humps’ seen in Seyfert spec- 
tra above ~7 keV and peaking near 20-30 keV. 

Another signature of reflection is the Fe Ka emission line 
commonly seen at 6.4 keV; this has proved to be a valu- 
able diagnostic of the accreting material near black holes. 
A few sources show evidence for a relativistically broad- 
ened component that is likely produced in the inner accre- 
tion disk, its profile sculpted by gravitational redshifting 
and relativistic Doppler effects (e.g., Fabian et al. 2002). 
However, a narrow Fe Ka component is a much more com- 
mon feature in Seyfert Is (e.g., O’Brien 2001, Yaqoob et 
al. 2001). FWHMs of several thousand km s _1 are typ- 
ical (e.g., Yaqoob & Padmanabhan 2004). The narrow 
component is generally thought to originate in Compton- 
thick material far from the central black hole, such as the 
outer accretion disk, the putative molecular torus invoked 


in standard unification models, or Compton-thick gas com- 
mensurate with the Broad-Line Region (BLR). 

At the same time, there is strong evidence from X-ray 
and UV grating observations for the presence of ionized 
material in the inner regions of a large fraction of AGN 
(e.g., Blustin et al. 2005). High-resolution spectroscopy 
shows the gas is usually outflowing from the nucleus; typi- 
cal velocities are ~ a few hundred km s _1 . Absorption due 
to a broad range of ionic species is commonly seen. For 
many sources, the relative line strengths argue for several 
different photoionized absorbing components, as opposed 
to a single absorber, along the line of sight. In the Fe K 
bandpass, absorption at 6.7 keV has been seen in NGC 
3783 with XMM-Newton (Reeves et al. 2004), consistent 
with absorption by Fe XXV. Additionally, absorption fea- 
tures near ~7-8 keV have been detected in PG and BAL 
quasars; such features have been attributed to strongly- 
blueshifted, high-ionization Fe K-shell absorption at near- 
relativistic (~0.2c) velocities (Reeves et al. 2003; Pounds 
et al. 2003a, 2003b; Chartas et al. 2002). These features 
may be signatures of high-velocity accretion disk winds. 

IC 4329a is a well-studied, nearby (z = 0.01605) Seyfert 
1.2 nucleus embedded in a nearly edge-on host galaxy. 
It has consistently remained one of the X-ray brightest 
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Seyfert galaxies for at least a decade. Using simultaneous 
ASCA and RXTE data, Done, Madejski & Zycki (2001) 
found the 6.4 keV Fe K a core to be moderately broad- 
ened (~20,000 km s -1 FWHM). The emission was param- 
eterized using a relativistic diskline model, with emission 
extending inward only to 30-100 gravitational radii R g . 
Using a 60 ksec CTiandro-HETGS observation performed 
in 2001, McKernan & Yaqoob (2004) detected not only a 
narrow 6.4 keV core, but also confirmed the presence of 
emission near 6.9 keV. This double-peaked complex could 
be fit by several competing models, included dual Gaus- 
sians, dual disklines, or a single diskline. Steenbrugge et 
al. (2005; hereafter S05) observed IC 4329a with XMM- 
Newton for 136 ksec in 2003 (see §2 for details). Their anal- 
ysis concentrated primarily on the soft X-ray RGS data, 
which revealed evidence for several absorbing components, 
including neutral absorption intrinsic to IC 4329a’s host 
galaxy, and a four-component warm absorber spanning a 
range of ionization states (the log of the ionization param- 
eter £ ranged from -1.4 to +2.7). Evidence for absorption 
due to local, z=0 hot gas was also detected. The EPIC-pn 
data (see below) also detected two emission peaks near 6.4 
and 6.9 keV (rest frame), identified as the Fe Ka core and 
a blend of Fe I K/3 and Fe XXVI emission, respectively. 
S05 reported that the 6.4 keV line was narrow and did 
not vary with time, consistent with an origin far from the 
central black hole. 

In this paper, we aim to not only parameterize the Fe 
K emission due to accreting material in IC 4329a, but also 
to search for absorption features in the Fe K bandpass 
that could be indicative of outflowing material. To this 
end, we have re-analyzed the XMM-Newton pn spectrum, 
concentrating on the Fe K bandpass. We also augment the 
XMM-Newton data with spectra derived from RXTE mon- 
itoring data to explore time-resolved spectroscopy over a 
wide range of time scales, search for variability in the Fe K 
core flux, and parameterize the relation between the core 
and the flux of the X-ray continuum. The XMM-Newton 
and RXTE observations and data reduction are described 
in §2. Spectral fits to the RXTE data are described in §3. 
Spectral fits to the Fe K bandpass of the XMM-Newton 
EPIC data are described in §4. Time-resolved spectral fit- 
ting to the RXTE and XMM-Newton data are described 
in §5. The results are discussed in §6, followed by a brief 
summary in §7. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. XMM-EPIC 

IC 4329a was observed by XMM-Newton during revolu- 
tion 670 on 2003 August 6-7, for a duration of 136 ksec. 
This paper uses data taken with the European Photon 
Imaging Camera (EPIC), which consists of one pn CCD 
back-illuminated array sensitive to 0.15-15 keV photons 
(Struder et al. 2001), and two MOS CCD front-illuminated 
arrays sensitive to 0.15-12 keV photons (MOSl and 
MOS2, Turner et al. 2001). Data from the pn were taken 
in Small Window Mode, data from the MOSl was in 
PrimePartialW3/large window mode, and data from the 
MOS2 was in PrimePartialW2/small window mode. The 
thin filters were used for the pn and MOS 2 cameras; the 
MOSl used the medium filter. Spectra were extracted us- 
ing XMM-SAS version 6.50, and using standard extraction 


procedures. For all three cameras, source data were ex- 
tracted from a circular region of radius 40"; backgrounds 
were extracted from circles of identical size, centered ~3' 
away from the core. Hot, flickering, or bad pixels were 
excluded. Data were selected using event patterns 0-4 for 
the pn. To reduce pile-up, only pattern 0 events were ex- 
tracted for the MOS cameras. However, we found (e.g., 
using the SAS task EPaplot) that the MOSl data about 
2 keV were severely piled-up; the level of pile-up ranged 
from ~5-10% at 2-3 keV to ~50% at 10 keV. We do not 
consider the MOSl data further in this paper. The MOS2 
data above 6 keV suffered from a small level of pile-up, 
ranging from ~3% at 6 keV to 10% at 10 keV. This reduc- 
tion yielded 133.3 ksec of pn data, starting at 2003 Aug 
06 at 06:57 UTC, and 132.9 ksec of MOS2 data, start- 
ing at 2003 Aug 06 at 06:51 UTC. After correcting for 
deadtime effects, the final exposure times for the pn and 
MOS2 were 93.2 ksec and 128.8 ksec, respectively. This 
paper will concentrate primarily on data obtained with 
the higher signal-to-noise pn, though MOS2 data were also 
analyzed for consistency. Background-subtracted pn light 
curves have already been presented in S05. We found a 
mean 2-10 keV pn count rate of 6.4 cs' 1 , corresponding 
to a flux of 8.9 x 10 -11 erg cm -2 s -1 . 

Since background flares due to proton flux tend to have 
hard spectra, the 10-12 keV band of the pn is the band 
most sensitive to them. The 10-12 keV pn background 
light curve, plotted in log-space, is shown in Figure 1. To 
test if these flares had any impact on the pn source spec- 
trum, we filtered out data collected during times when 
the 10-12 keV background light curve exceeded a count 
rate of B + 2a (equal to 0.3 c s _1 ), where B is the mean 
background rate and a is the standard deviation of the 
light curve. This screening yielded a new exposure time 
of 71.6 ksec. We found that the only major difference 
was that the screened spectra tended to be fit with a pho- 
ton index approximately 0.1 steeper than the unscreened 
spectra, but all other parameters remained virtually un- 
changed. In particular, the Fe K profile, the focus of this 
paper, remained unchanged. We therefore concentrated on 
the unscreened data in all analysis below to increase the 
signal-to-noise, though we checked the screened data for 
consistency. However, we note that the brightness of the 
sources could have an impact on the degree of difference 
between the screened and unscreened spectra; IC 4329a is 
an X-ray bright AGN, and for fainter sources, such screen- 
ing could potentially have a much larger impact than for 
this dataset. We therefore recommend that the commu- 
nity explore this type of screening as an additional precau- 
tion against source contamination from background flares, 
especially for fainter sources. 

2.2. RXTE observations and data reduction 

IC 4329a has been monitored by RXTE once every 4.3 d 
since 2003 Apr 8 (hereafter referred to as long-term sam- 
pling). In this paper, we include data taken up until 2005 
Oct 02. Due to sun- angle constraints, no data were taken 
from 2003 Oct 03 - 2003 Nov 20 or from 2004 Oct 04 - 2004 
Nov 21. IC 4329a was subjected to intensive RXTE moni- 
toring, once every third orbit (17 ksec) from 2003 Jul 10 to 
2003 Aug 13 (hereafter referred to as medium-term sam- 
pling; ’’short-term” sampling refers to the XMM-Newton 
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long-look detailed above). Each RXTE visit lasted ~1 
ksec. 

2.2.1. PCA data reduction 

Data were taken using RXTE 1 s proportional counter ar- 
ray (PCA; Swank 1998), which consists of five identical col- 
limated proportional counter units (PCUs). For simplic- 
ity, data were collected only from PCU2 (PCUO suffered 
a propane layer loss in 2000 May; the other three PCUs 
tend to suffer from repeated breakdown during on-source 
time). Count rates quoted in this paper are per 1 PCU. 
The data were reduced using standard extraction meth- 
ods and FTOOLS V5.3.1 software. Data were rejected if 
they were gathered less than 10° from the Earth’s limb, 
if they were obtained within 30 min after the satellite’s 
passage through the South Atlantic Anomaly, if ELEC- 
TRON > 0.1, or if the satellite’s pointing offset was 
greater than 0?02. 

As the PCA has no simultaneous background monitor- 
ing capability, background data were estimated by using 
PCABACKEST v3.0 to generate model files based on the 
particle-induced background, SAA activity, and the dif- 
fuse X-ray background. The ’L7-240’ background models 
appropriate for faint sources were used. This background 
subtraction is the dominant source of systematic error in 
RXTE AGN monitoring data (e.g., Edelson & Nandra 
1999). Counts were extracted only from the topmost PCU 
layer to maximize the signal-to-noise ratio. Spectra from 
individual visits were merged to form a long-term time- 
averaged spectrum (all data) and a medium-term time- 
averaged spectrum (intensive monitoring only); total ex- 
posures were 278.1 ksec and 100.4 ksec, respectively. 

Response matrices were generated using PCARSP v.8.0 
over several periods between 2003 and 2005. The PCA re- 
sponse evolves slowly; due to the gradual leak of propane 
into the xenon layers, the spectral response hardens slowly 
with time. For example, application of a 2003 matrix to 
IC 4329a data observed in 2005 (or vice versa) results in a 
shift in photon index of only ~0.04. For long-term time- 
averaged spectral analysis below, a 2004 response matrix 
was used. 

2.2.2. HEXTE data reduction 

The High-Energy X-Ray Timing Experiment (HEXTE) 
aboard RXTE consists of two independent clusters (A and 
B), each containing four scintillation counters (see Roth- 
schild et al. 1998) which share a common 1? FWHM field 
of view. Source and background spectra were extracted 
from each individual RXTE visit using Science Event data 
and standard extraction procedures. The same good time 
intervals used for the PCA data (e.g., including Earth el- 
evation and SAA passage screening) were applied to the 
HEXTE data. To measure real-time background measure- 
ments, the two HEXTE clusters each undergo two-sided 
rocking to offset positions, in this case, to 1.5° off-source, 
switching every 32 seconds. There is a galaxy cluster 
(Abell 3571, z = 0.039) located about 2° south of IC 4329a, 
at R.A. = 13h 47.5m, deck = -32° 52 m. This source is de- 
tected in the RXTE all-sky slew survey (XSS; Revnivtsev 


et al. 2004), which shows the 8-20 keV flux of this source 
to be about half that of IC 4329a. This suggests that the 
presence of A3571 could affect HEXTE background deter- 
mination for any HEXTE cluster which rocks to within 
roughly 1-2° of the center of A3571 (accounting for the 
extended size of the galaxy cluster and the ~1° FWHM 
collimator response of the HEXTE clusters). However, 
BeppoSAX- PDS observations have shown no detection of 
A3571 above 15 keV (Nevalainen et al. 2004). The 8-20 
keV emission seen by the XSS must therefore be emission 
only between 8 and <15 keV, and so the presence of A3571 
can be safely ignored as far as contaminating the HEXTE 
background data is concerned. Cluster A data taken be- 
tween 2004 Dec 13 and 2005 Jan 14 were excluded, as the 
cluster did not rock on/off source. Detector 2 aboard Clus- 
ter B lost spectral capabilities in 1996; these data were ex- 
cluded from spectral analysis. Cluster A and B data were 
extracted separately and not combined, as their response 
matrices differ slightly. Deadtime corrections were applied, 
and the individual spectra were merged for each cluster to 
form a long- and medium-term time-average spectrum in 
the same manner as the PCA data. The net source ex- 
posure times for the long-term spectrum were 100.5 ksec 
(cluster A) and 101.8 ksec (cluster B). Net source expo- 
sure times for the medium-term spectrum were 37.0 ksec 
(A) and 36.5 ksec (B). 

3. pca/hexte time-average spectral analysis 

We first discuss spectral fits to the long- and medium- 
term RXTE PCA and HEXTE data. The energy resolu- 
tion of the PCA is low, but we discuss these data first for 
the purpose of characterizing the broadband hard X-ray 
continuum and constraining the strength of the Compton 
reflection component, and applying those constraints to 
the the high-resolution XMM-Newton EPIC data. Addi- 
tionally, model fits to the long- and medium-term time- 
averaged spectra were used as templates when analyz- 
ing the time-resolved spectra, as described in §5. The 
medium-term data constitutes a large fraction (36%) of 
the total long-term data, so the two spectra are not inde- 
pendent. 

PCA data below 3.5 keV were discarded in order to 
disregard PCA calibration uncertainties below this en- 
ergy and to reduce the effects of the warm absorber (e.g., 
KS05). For the long-term data, PCA data above 38 keV 
were excluded as the source counts become dominated by 
statistical uncertainties at higher energies; in the medium- 
term data, PCA data above 30 keV were excluded. All 
PCA data from 24-32 keV were also ignored due to the 
presence of a feature in the background that is not prop- 
erly modeled by the current PCA calibration for faint 
sources (C. Markwardt, priv. comm.). 1 HEXTE data be- 
low 25 keV were excluded, as the responses of clusters A 
and B diverge somewhat for faint sources (N. Shaposh- 
nikov, priv. comm.); there is good agreement above 25 
keV. HEXTE data above 100 keV were excluded as the 
source count rate becomes dominated by statistical uncer- 
tainties at higher energies. For the long-term time-average 
spectrum, HEXTE data above 50 keV were grouped as fol- 
lows: channels 51-60, channels 61-75, and channels 76 and 


1 As of the time this paper was submitted, the PCA group was preparing to release new background models; testing on IC 4329a showed that 
the 24-32 keV background features are removed; these data will be included if the new backgrounds are released to the public before the paper 
is accepted. 
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above were binned in groups of 2, 3, and 4, respectively. 
For the medium-term data, data above 39 keV were also 
binned: channels 39-50 and channels 51 and above were 
binned in groups of 3 and 8, respectively. 

All spectral fitting for both RXTE and XMM-Newton 
data was done with XSPEC v. 11. 3. 21 (Arnaud 1996). Er- 
rors quoted for spectral fit parameters are 90% confidence. 
All line energies quoted are for the rest frame unless oth- 
erwise indicated. All RXTE and XMM-Newton spectral 
fits below included a WABS neutral absorption component. 
The Galactic column towards IC 4329a is 4.42 xlO 20 cm -2 
(Dickey & Lockman 1990), but previous studies have sug- 
gested the presence of additional neutral absorption (e.g., 
Gondoin et al. 2001). As we did not analyze <2.5 keV 
XMM-Newton pn data or <3.5 keV RXTE PCA data, this 
analysis is relatively insensitive to the precise value of the 
total (Galactic + intrinsic) neutral column density. We 
thus left the neutral column density as a free parameter in 
all fits below. 

We first discuss joint fits to the long-term PCA/HEXTE 
spectra, plotted in Figure 2. In these fits, a constant coeffi- 
cient was included to account for minor normalization off- 
sets between the PCA and HEXTE. The first model tested, 
a simple power-law modified by cold absorption (Model 1), 
yielded large residuals near 6 keV due to the presence of 
the Fe Ka line; x 2 was 2129.17/128 degrees of freedom 
(. dof ). The addition of a Gaussian emission component 
(Model 2) at 6.38 keV, the Fe Ka line, proved significant at 
greater than 99.99% confidence in an F-test, since x 1 /dof 
fell to 586.81/125. However, there still remained broad 
residuals due to the presence of the Compton reflection 
component. Model 3 accounted for the presence of this 
component by using a PEXRAV component to model reflec- 
tion of an underlying power-law component off optically 
thick and neutral material (Magdziarz & Zdziarski 1995). 
The inclination angle i was initially frozen at 30°, and the 
power-law cutoff energy was kept fixed at 270 keV as per 
the results of Perola et al. (1999). Elemental abundances 
were kept fixed at solar. Adding the PEXRAV component 
proved significant at greater than 99.99% confidence in 
an F-test, since x 2 /dof fell to 208.14/123. The photon 
index T was 1.894io.oii- The reflection fraction R (ra- 
tio of the reflection normalization at 1 keV to that of the 
power law) was 0.51 ±0.04. (With HEXTE data omitted, 
R from PCA data alone is slightly higher at 0.55±0.07; 
T = 1.908 ±0.015.) For an near face-on inclination (cos(i) 
= 0.95), x 2 increased by 0.1, and R dropped slightly to 
0.48±0.04. With the constant coefficient of the PCA fixed 
at 1.0, the constants for the HEXTE A and B spectra were 
0.90 and 0.92, respectively. Table 1 summarizes the model 
parameters (including line emission parameters) for Model 
3. Figure 2 shows the data/model residuals for the three 
models discussed. 

Given the limit of PCA calibration (data/model resid- 
uals ~2% 2 ) and the high signal-to-noise ratio of this 
spectrum, systematic errors associated with background 
modeling and the PCA’s spectral response dominated 
the uncertainties. Finding a statistically-acceptable time- 
averaged fit (with reduced x 2 ~ 1) was thus unlikely even 
if the model was representative of the intrinsic spectrum of 
the source. Additionally, given the low energy resolution 


of the PCA (~1 keV at 6 keV), the PCA is largely insensi- 
tive to the detailed shape of the Fe Ka line; there were no 
other obvious residuals seen in the Fe K bandpass. For ex- 
ample, Fe K/3 emission is clearly detected in the EPIC-pn 
spectrum at 6.90 keV as discussed below, but this is not 
obvious from the PCA residuals; adding a second Gaussian 
with rest-frame energy centered at 7.06 keV was only sig- 
nificant at 90% confidence. Additionally, using a diskline 
model the emission profile (as will be discussed for the 
EPIC-pn data below) did not provide any improvement in 
the fit. Model 3 is therefore adopted as a satisfactory fit 
to the long-term time-averaged PCA/HEXTE spectrum. 
Since the primary goal of the PCA/HEXTE fitting was to 
constrain the photon index T and the reflection strength R, 
we will study the more detailed models (adding a second 
Gaussian to model the Fe K/3 emission; modeling Compton 
reflection off ionized material; modeling the Fe K emission 
using relativistic disklines) only using the high resolution 
EPIC-pn data (§4), checking the PCA spectrum for con- 
sistency. 

Spectral fitting to the medium-term time average 
PCA/HEXTE spectrum proceeded in a manner identical 
to the long-term spectrum, testing the same three models, 
and with very similar results, although x 2 /dof was lower 
due to lower signal-to-noise. The results of the Model 3 
fit are summarized in Table 1. Compared to the long- 
term time average spectrum, the medium-term T and R 
are both found to be slightly lower. However, we note the 
similarity in line emission parameters. Furthermore, the 
fact that the same type of model is representative of both 
spectra is consistent with the notion that the general form 
of the broadband X-ray emission did not change over time. 

4. PN FITS 

We now focus on spectral fits to the time-average EPIC 
pn spectrum, shown in Figure 3. The best-fitting model 
will be used as a template for time-resolved spectroscopic 
studies in §5. We ignore <2.5 keV data and concen- 
trate only on hard X-ray emission. Additionally, the ef- 
fective area of the pn drops off rapidly below 2.5 keV. All 
model fits included absorption by neutral material (WABS 
model), with the column density left as a free parame- 
ter. The unscreened data was used in order to maximize 
the total good exposure time. §4.1 describes preliminary 
fits to the data, including attempts to model the Fe K 
emission using simple Gaussians. §4.2 describes fits using 
relativistic diskline models. In §4.3, we discuss a narrow 
absorption feature seen at 7.68 keV, and discuss its detec- 
tion significance as well as its possible origin. In §4.4, we 
show that the pn spectrum is consistent with the MOS2 
and simultaneous RXTE spectra. 

4.1. Fe Ka and Kf3 emission: preliminary fits 

Model 1 , consisting of a simple power law, yielded a sta- 
tistically unacceptable fit with large data/model residuals 
which clearly indicated the presence of an emission line 
near 6.4 keV (the Fe Ka core) and a rise in the continuum 
above 7 keV (Compton reflection hump). Data/model 
residuals for the best-fit model are shown in Figure 3. For 
Model 2, a narrow Gaussian was added to model the Fe Ka 
line; the best-fitting energy is 6.40±0.01 keV, and the core 


2 See http://lheawww.gsfc.nasa.gov/users/keith/rossi200Q/energy_response.ps 
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width <ti is 112+19_i5 eV. However, there were still sig- 
nificant data/model residuals near 7 keV suggesting addi- 
tional line emission, and residuals in the continuum above 
7 keV which suggested that the Compton reflection hump 
is detected independently of the Fe Ka core. For Model 3, 
a PEXRAV component was added, with R fixed at 0.51, the 
value from the long-term RXTE PCA/HEXTE fits above. 
The high-energy cutoff was fixed at 270 keV (e.g., Perola 
et al. 1999); solar metal abundances and an inclination 
of 30° were assumed. The best-fitting model had x 2 /do f 
= 1865.73/1893, but data/model residuals, shown in Fig- 
ure 4, indicated strong, narrow excess emission near 7.0 
keV as well as a deficit near 7.7 keV which suggested an 
absorption feature of some kind. The best-fit spectral pa- 
rameters for Models 1-3 are listed in Table 2. As an inde- 
pendent confirmation of the edge at 7.1 keV, we fitted the 
pn data with an edge instead of a PEXRAV component. We 
found an edge energy of 7.16+0.07 and an optical depth 
of 0.05+0.01, similar to the results found by Gondoin et 
al. (2001). 

We hypothesized that the emission near 7.0 keV (the 
“blue peak”; “red peak” refers to the 6.4 keV line) could 
be due to Fe K/3, H-like Fe K, or a blend of both. Accord- 
ingly, Model 4 featured the inclusion of a second Gaussian 
component near 7.0 keV. The best fit had a blue peak 
with energy centroid at 6.931 q °q keV and a width <r 2 of 
121±I! eV. The energy centroid of the Ka Gaussian did 
not change, though the Ka width oq decreased slightly to 
91+ 13 eV. The best-fit x 2 /dof was 1844.88/1890. Using 
an F-test to assess the significance of adding this compo- 
nent yielded a null hypothesis probability of 9.36xl0~ 5 . 
Data/model residuals are plotted in Figure 4; the best-fit 
spectral parameters are listed in Table 2. 

In this model, the red peak centroid energy, 6.39+0.01 
keV, is consistent with Ka emission from neutral (or at 
best mildly-ionized) Fe. The core width cq , 91+13 eV, cor- 
responds to a FWHM width of 210 ± 30 eV, or a FWHM 
velocity of 9700 ± 1400 km s -1 . This value is consistent 
with that derived using Chandra-HETGS, 151001 jqooo km 
s" 1 (Yaqoob & Padmanahban 2004). We note that the 6.4 
keV line is resolved at the pn resolution of ~140 eV: fixing 
<ti at 1 eV caused large data/model residuals to appear 
between 6. 0-6. 6 keV while x 2 increased by 91. Similarly, 
the blue peak is also resolved by the pn. Fixing <r 2 to 1 
eV caused x 2 /dof to increase to 1851.40/1891, significant 
at 99.0% confidence in an F-test. 

Given the energy resolution of the pn, it was not ob- 
vious from these parameters whether the emission line is 
due to neutral (or only mildly-ionized) Fe K/3 at 7.056 
keV, Fe XXVI at 6.966 keV, or a blend of both. For ex- 
ample, fixing cr 2 to be equal to oq, as would be expected 
if the 7.0 keV emission is due entirely to Fe K/3, resulted 
in very little change in x 2 (-0.86), and virtually no change 
in any other model parameter; i.e., it was not significant 
to thaw <j 2 . To test if the blue peak was consistent with 
being due solely to neutral Fe K/3 emission, we adjusted 
Model 4 as follows: we fixed the blue peak’s centroid en- 
ergy at 7.06 keV and fixed the ratio of the K/3/Ka fluxes 
at 17/150. The width c 2 was kept tied to that of the Fe 
Ka core, oq. The best fit model (Model 5) showed vir- 
tually no change in the 6.4 keV core or continuum model 
parameters compared to Model 4 (see Table 2). However, 


the fit was slightly worse: x 2 /dof was 1852.78/1893, and, 
as shown in Figure 4, small data/model residuals appeared 
near 6.8 keV, signifying excess unmodeled emission on the 
red side of the blue peak. 

As shown in Figure 4, the data/model residuals for Mod- 
els 4 and 5 suggested unmodeled emission on the red side of 
the Ka line 6.1 keV, prompting us to test for the presence 
of a Compton-scattering shoulder. Such a feature might 
be expected at 6.24 keV due to Fe Ka photons Compton- 
scattering on electrons in the medium in which the line 
is formed. The ratio of the shoulder and core intensities 
indicates scattering in, and an origin for the Fe Ka line 
in, Compton-thick or Compton-thin material (ratio above 
or below ~0. 1-0.2, respectively; see Matt 2002 for further 
details). Building on Model 4, the better-fitting of the 
double-Gaussian models, we added a third Gaussian emis- 
sion component with centroid energy fixed at 6.24 keV and 
a core width tied to that for the 6.4 keV core (Model 6A). 
The best-fit model had a x 2 /dof of 1839.03/1889, and an 
equivalent width EW of 11+8 eV. An F-test showed it was 
significant at 98.6% to add this component, however, the 
change in x 2 compared to Model 4 was only 5.85. Thawing 
the centroid energy (Model 6B) resulted in a modest im- 
provement in the fit; x 2 /dof dropped to 1833.68/1888 as 
the centroid energy moved to now 6. 131^25 keV. Best-fit 
spectral parameters for Model 6A and 6B are listed in Ta- 
ble 3. According to an F-test, it was significant at 98.1% 
confidence to thaw the centroid energy. Considering the 
pn energy resolution, this feature does not seem entirely 
consistent with emission from a Compton shoulder. 

4.2. Fe Ka and K(3 emission: diskline fits 

We now discuss whether relativistic diskline models 
(Laor 1991) can describe the Fe K profile of IC 4329a. 
Done et al. (2000) noted from the ASCA observation that 
IC 4329a does seem to have a slightly broadened red wing, 
though not as broad as that of the “archetypal” broad line 
source, MCG-6-30-15 (Tanaka et al. 1995); see Figure 2 of 
Done et al. (2000). Their best-fitting models for the broad 
emission suggested that there was not significant line emis- 
sion at radii smaller than ~30-100 R g (R g = GMbh/c 2 ). 
However, compared to the ASCA data, the pn spectrum 
reveals additional emission and absorption features in the 
Fe K band, prompting us to probe if diskline models can 
successfully account for all of the emission/absorption fea- 
tures in the Fe K bandpass. 

As seen by McKernan & Yaqoob (2004), the approxi- 
mate heights above the power-law continuum of the red 
and blue peaks in the 60 ksec HEG spectrum were ap- 
proximately equal, with the blue peak possibly slightly 
higher. This prompted McKernan & Yaqoob (2004) to 
test whether the two peaks could be horns from the same 
diskline model. Their best-fit model was one with an incli- 
nation angle of 24j/®°, rest-frame line energy of 6.74^° ^ 
keV, a flat radial emissivity profile (/3 < 0.7, where the 
radial emissivity per unit area is quantified as a power 
law, r _/? ), and emission between ~6 and 70 R g . However, 
the current EPIC-pn spectrum clearly shows that the EW 
and height above the power-law continuum of the red peak 
are much greater than those of the blue peak. A rela- 
tivistic diskline model with near face-on inclination, e.g., 
<10°, can typically produce only one peak, and thus one 


6 


Fe K Emission and Absorption in IC 4329a 


diskline cannot account for both peaks simultaneously. A 
diskline model with intermediate or high inclination angles 
is required to produce two peaks; usually the blue horn is 
much larger than the red horn unless the disk is truncated 
at ~tens of R g and/or the radial emissivity profile is rela- 
tively flat. We attempted to fit such a model to the Fe K 
profile in the pn data. Model 7 consisted of a power-law, a 
PEXRAV component, and a LaoR diskline component. We 
tested inclination angles of 30°, 50°, and 70°. The outer 
radius R ou t was kept fixed at 400 R g . The inner radius 
R m was kept free, but in all fits, R in went to very large 
values, usually 250-350 R g , with very poor constraints. 
Several values of the emissivity index /3 were tested (1,2, 
3, 5, 7), but all fits yielded poor or no constraints on /3. 
In all cases, there were large systematic data/model resid- 
uals across the Fe K bandpass, as shown in Figure 4 (see 
also Table 4). The values of x 2 for best-fit models were 
typically 1920-2140 for 1893 dof, significantly worse fits 
than e.g., the double-Gaussian model. We conclude that 
a single-diskline model is an inaccurate description of the 
data. 

The unmodeled excess emission near 6.1 keV and 6.8 
keV in the double-Gaussian model (Model 5) potentially 
signified that the Ka and K (3 emission components could 
each be independently modeled better using a redward 
asymmetric, moderately relativistic diskline profile instead 
of a symmetric Gaussian profile. McKernan & Yaqoob 
(2004) tested dual-diskline models for the HEG spectrum; 
the best-fit model was one with a nearly face-on inclination 
i < 6°, a moderately steep emissivity index (3 — 2.4±0.2, 
and emission spanning the radii between 6 R e (fixed) to 
600 R g . S05’s best fit double-diskline model was found to 
be consistent with this; (3 was < 1.4, an inclination angle 
<17°, a profile that is very similar to the double-Gaussian 
model. We also fit a double-diskline model to the data, 
with i, /3, R[n , and R out for the red and blue disklines 
set equal to each other, and the normalization of the blue 
diskline equal to 17/150 that of the red diskline. R ou t was 
kept fixed at 400E 9 ; it was not significant to thaw this pa- 
rameter. Our best fit model (Model 8, see Table 5) agrees 
generally well with that of S05; we found i = ll J 1 ®, (3 = 
2.1T0.3, and R m — 13t (1,5 R g ■ The energy of the red peak 
was 6.441 o qj keV, more consistent with Ka emission from 
mildly ionized Fe than from neutral Fe. The blue peak en- 
ergy was 6.98±0.09 keV, consistent with K/3 emission from 
neutral or mildly-ionized Fe. The equivalent width EW^ a 
of the Fe Ka line in this model is 86*g 7 eV. The values 
of EWK a and Compton reflection strength R were consis- 
tent with the prediction of George & Fabian (1991) that, 
for reflection off neutral material, the reflection strength 
R will be equal to EW^ a / 150 eV. 

As shown in Figure 4, the dual-diskline model accounts 
well for all of the excess emission- like residuals near 6.0 and 
6.8 keV, but it does not account for the 7.7 keV absorption- 
like residuals, which still persist. The value of x 2 /dof for 
this model was 1835.02/1890, i.e. , the value of reduced 
X 2 is slightly lower than that for the Compton shoulder 
model, 6B. The fact that the 6.1 and 6.8 keV residuals 
are now well modeled suggests that this model is a better 
description of the data compared to the double Gaussian 
models. To quantify this statement, we replaced only the 
red peak Gaussian in Model 5 with a Laor component and 


refitted, x 2 /dof fell to 1843.92/1890, significant at 97.2% 
confidence in an F-test compared to Model 5. Replacing 
only the blue peak Gaussian with a Laor component and 
refitting yielded x 2 /do/=1846. 48/1889, significant at only 
83.1% confidence in an F-test compared to Model 5. Ad- 
ditionally, the fact that Model 8 takes care of the ~6.0 
keV residuals suggests it is a better fit than models with a 
Compton shoulder. In fact, adding a Gaussian in emission 
at 6.24 keV to Model 8, we found only an upper limit of 6 
eV to Compton shoulder emission, e.g., we conclude that 
do not detect the Compton shoulder with high significance. 

A Laor profile assumes a maximally-rotating Kerr 
black hole. Substituting the Laor model components 
with DISKLINE model components and assuming a non- 
rotating Schwarzschild black hole (Fabian et al. 1989), into 
the model yields a very similar fit to Model 8 (x 2 /dof = 
1839.36/1890); one would not expect significant differences 
in the line profile unless there was significant emission 
originating from within a few R g of the black hole. Al- 
most all DISKLINE best-fit parameters were virtually iden- 
tical to the Laor parameters, though the inclination was 
slightly higher ( i = 14^®°). For the remainder of this pa- 
per, “diskline” will denote the Laor model only. 

Another possible diskline model is one wherein the 
Doppler horns are very close together and unresolved by 
the pn. A flat emissivity profile and an inclination angle 
of ~20-30° can yield such a profile for certain ranges of 
/ 3 and R in . Additionally, the blue horn in this profile is 
slightly higher than the red horn; the fact that the best- 
fitting rest-frame energy for the Fe Ka diskline in Model 
8 was 6.44±0.01 suggests that this model is worth explor- 
ing. Inclination angles much higher than ~30° yield line 
profiles with Doppler peaks that could be resolved by the 
pn. We tried applying such a model to the data, fixing the 
inclination angle to values in the range 20-30 °. However, 
these models could not account for the peak of the emission 
and the broad red wing simultaneously. Values of x 2 /dof 
generally stayed about 1900/1891, and we conclude that 
this model does not describe the emission profile as well 
as Model 8 above. 

Finally, we explored the possibility that there could be 
additional contributions to the emission profile from nar- 
row Fe Ka and K/3 lines, as would be expected if there 
is non-negligible Fe K emission originating in material far 
from the black hole. We added two narrow Gaussians to 
Model 8 at the rest-frame energies for Fe Ka and K/3, with 
the K/3 normalization set to 17/150 times that of the Ka 
line. We first considered the possibility that each broad 
component is described by a low-inclination diskline com- 
ponent, similar to that used in Model 8. We kept the 
widths of both peaks fixed at 50 eV. The best-fit model 
(referred to as Model 9A hereafter) is one with diskline 
parameters nearly identical to those in Model 8, with the 
Ka and K/3 lines contributing only a small amount, and 
formally consistent with upper limits only in EW. The 
best-fit parameters are listed in Table 6. Compared to 
Model 8, x 2 decreases by only 1.2 for one additional dof. 
Data/model residuals look identical to those in Model 8 
and are not plotted. The maximum allowable EW for 
the narrow Ka line, without significantly changing any of 
the diskline parameters, is 27 eV; for this case, the EW 
of the broad Ka line is 54 eV. We conclude that if the 
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two broad components are described by low-inclination, 
moderately-relativistic disklines, then the bulk or entirety 
of the emission must be from the broad component, but a 
contribution from two narrow components cannot be ruled 
out. 

We next considered the possibility that for each peak, 
a high-inclination broad line component is present and 
responsible primarily for the red-wing emission, with a 
narrow component accounting for the bulk of the peak 
emission very close to 6.40 or 7.06 keV. We kept the rest- 
frame energies of the Fe p broad and narrow components 
fixed at 7.06 keV. Parameters for the best-fit model (hence- 
forth referred to as Model 9B) are listed in Table 6. The 
best-fitting x was nearly identical to that for Model 8. 
Data/model residuals look virtually identical to those in 
Model 8 and are not plotted. The width <r of the narrow 
Ka line, 66+22 e ^> corresponds to a FWHM velocity of 
7000 km s -1 . For Keplerian motion, this corresponds to a 
radius of ~900f?sch (1 #ScH=2GMbh/c 2 )- 

Given that models with relativistic disklines (models 8, 
9A, and 9B) provide a superior fit to the emission pro- 
file compared to the double-Gaussian models (4 and 5) 
or models with double-Gaussians plus a Compton shoul- 
der, particularly in accounting for the presence of a Ka 
red wing, we can conclude that some sort of diskline is 
indeed required to best model the emission. However, a 
degeneracy in modeling the profile is apparent: models 
containing two broad lines and two narrow lines can be fit 
the observed profile of IC 4329a equally well using either a 
low-inclination diskline (plus little or no narrow emission) 
or by narrow components plus high-inclination disklines. 

To further quantify the presence of emission features 
in the Fe K bandpass in addition to the Fe Ka and K p 
profiles, we added a narrow Gaussian to Model 8, and, us- 
ing the STEPPAR command in XSPEC, derived confidence 
contours of line intensity versus rest-frame energy. We 
stepped through values of the energy centroid in the range 
5.60-8.30 keV in increments of 0.05 keV and flux range 
0-3 x 10~ 5 ph cm -2 s -1 in increments of 1 x 10 -6 ph cm -2 
s -1 . The width cr was kept fixed at 100 eV. The results are 
shown in the top panel of Figure 5, and demonstrate that 
there were no additional obvious signatures of emission, 
including those from the Fe XXV triplet or from Fe XXVI 
Lya at their rest-frame energies. Similarly, adding an in- 
verted Gaussian to quantify absorption, searching over the 
flux range 0 to -2.2xl0~ 5 ph cm -2 s -1 showed that there 
was no strong evidence from these plots for strong reso- 
nance absorption at the rest-frame energies of Fe XXV or 
Fe XXVI (see Figure 5, bottom panel). The only strong 
narrow absorption feature was the 7.7 keV feature. 

To further additional quantitative limits on K-shell 
Fe XXV and Fe XXVI emission and resonance absorp- 
tion, we modified Model 8 by adding a narrow (a = 65 
eV) Gaussian (inverted for absorption). The energy cen- 
troid was fixed at 6.70 or 6.97 keV for He- or H-like Fe; 
though given the energy resolution of the pn, it would not 
have been possible to determine whether it was the for- 
bidden, intercombination, or resonance line that was re- 
sponsible for the He-like emission. In all cases, we found 
that adding a Gaussian increased y 2 , and found only up- 
per limits to emission or absorption. The upper limits 
on the absolute values of the intensity ( EW ) for Fe XXV 


emission, Fe XXVI Lya emission, Fe XXV absorption, or 
Fe XXVI Lya absorption were 8.5x1 0~ 6 (8 eV), 13.2 x If) -6 
(15 eV), 2.4xl0~ 6 (3 eV), and 1.8xl(r 6 ph cm" 2 s' 1 (2 
eV), respectively, consistent with the results of Bianchi et 
al. (2005) for the 10 ksec XMM-Newton observation of IC 
4329a in 2001. 

From our analysis, we conclude that the sets of models 
containing dual-disklines (8, 9A, or 9B) are equally ade- 
quate in fully describing the Fe K emission profile of IC 
4329a, since there is no formal requirement for additional 
emission lines. 

4.3. The 7.68 keV narrow absorption feature 

We now turn our attention to modeling the narrow ab- 
sorption feature which appears at ~7.7 keV. This feature 
cannot be attributed to relativistic diskline emission, and 
appears in the data/model residuals independently of how 
the Fe K emission lines are parameterized. We assume that 
it has a physical origin distinct from that of the emission 
lines. §4.3.1 discusses simple model fits to the absorp- 
tion line. In §4.3.2, we show that this feature is likely not 
an artifact of the background, and in §4.3.3, we discuss 
the statistical significance of detecting this feature, and 
present the results of a Monte Carlo simulation analysis 
which show that the feature is almost certainly not due 
to photon noise. Based on these results, we can safely as- 
sume the feature is intrinsic to the source, and in §4.3.4, 
we discuss possible origins for the line. 

4.3.1. Gaussian model fit for the line 

We started with Model 8, the power-law + PEXRAV + 
double-diskline model (though the reader is reminded that 
Models 9A and 9B provide an equally acceptable fit to the 
emission profile, and the spectral-fit results for the absorp- 
tion line are independent of which model is used to de- 
scribe the emission profile). We first tested the hypothesis 
that the feature was indicative of reflection off ionized ma- 
terial. First, we substituted the PEXRAV component with 
a PEXRIV component (Model 10); an increase in the ion- 
ization level of the reflecting material pushes the energy of 
the Fe K absorption edge from 7.1 keV up towards higher 
energies. However, using a PEXRIV component, with disk 
temperature T^sk kept fixed at 30,000 K, resulted in a best 
fit virtually identical to Model 8. The value of the ioniza- 
tion parameter £ (defined as AnF ion /n, with F ion defined as 
the 0.5-20 keV ionizing continuum flux and n is the density 
of the reflecting material) was 0 erg cm s -1 (upper limit of 
0.04 erg cm s -1 ). x 2 /dof increased to 1838.46/1889 with 
no improvement in the data/model residuals. Forcing the 
Fe K edge energy to be near 7A-7.7 keV by constraining £ 
to be near ~3-10 erg cm s -1 did not successfully remove 
the residuals. Thawing Tdisk provided no improvement to 
the fits. Models simultaneously incorporating PEXRAV and 
PEXRIV components, with absorption edges near 7.0 keV 
and 7.4-7.5 keV, respectively, also did not account for the 
7.7 keV residuals. In all cases, an Fe K absorption edge 
near 7.1 keV, consistent with reflection off of neutral ma- 
terial, was still required by the data (see, e.g., Gondoin et 
al. 2001). 

We tested the notion that the feature could be an ab- 
sorption edge, independent from Compton reflection. We 
returned to Model 8 (with PEXRAV), and included an edge 
component (Model 11). The best-fit model was one with 
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an edge energy of 7.17^0$ keV and an optical depth t 
of 0.03±0.01. However, the data/model residuals near 7.6 
keV were not affected. Attempts to force the edge energy 
closer to 7. 4-7. 6 keV were also unsuccessful in removing 
the data/model residuals. Spectral parameters for the best 
fit Models 10 and 11 are listed in Table 7. From the re- 
sults of the Model 10 and 11 fits, one can conclude that 
the feature is too narrow to be an absorption edge. 

We then attempted to model the feature using a sim- 
ple narrow, inverted Gaussian. We returned to Model 8, 
and added an inverted Gaussian near 7.7 keV (Model 12). 
The best-fit model successfully removed the data/model 
residuals, as shown in Figure 4. Spectral parameters for 
the best fit Model 12 are listed in Table 7. The best-fit 
rest-frame energy centroid was 7.68 tg ; ^ keV (correspond- 
ing to 7.56+0 03 keV in the observed frame). The best fit 
width a was 35+35 eV. The absolute value of the intensity 
was 9.6±3.4 x 10 -6 ph cm -2 s -1 for an \EW\ of 13±5 
eV. The value of x 2 /d°f was 1813.50/1887, the best-fit 
model so far. Compared to Model 8, x 2 dropped by 21.73, 
and an F-test for adding this component yielded F=7.54, 
with a null hypothesis probability of 5.2 xlO -5 . However, 
as discussed in §4.2.3, there was no a priori reason to ex- 
pect an absorption feature in this energy bandpass, and 
the F-test, when used in this “standard” manner, tends to 
overestimate the significance of detecting a feature. The 
statistical significance of detecting this component will be 
addressed in more detail in §4.2.3. 

Finally, having properly quantified the emission and 
absorption profile, we investigated the Fe abundance by 
studying the height of the 7.1 keV edge relative to the 
Compton reflection strength. In the PEXRAV model, the 
Fe abundance had been kept frozen at 1.0, using the An- 
ders & Grevesse (1989) tables. Keeping R frozen at 0.51 
and thawing the Fe abundance made no change in the fit 
statistic; the best-fit Fe abundance was l.OlJ ®. Using the 
Grevesse & Sauval (1998) tables, the best-fit Fe abundance 
was 1.5±0.7. 

4.3.2. PN background/instrumental effects 

The pn background spectrum shows no obvious features 
around observed-frame energies of 7.5-7. 6 keV. However, 
two instrumental features, a Ni Ka emission line at 7.48 
keV and a Cu Ka emission line at 8.05 keV, are expected 
(Katayama et al. 2004; Freyberg et al. 2004). These lines 
originate from the electronics board located below the 
CCD camera. However, there are no obvious features in pn 
background spectrum for IC 4329a at these energies (see 
Figure 6). One might initially suspect that the 7.56 keV 
(observed-frame) absorption feature could be an artifact 
resulting from the 7.48 keV instrumental line, but there 
are several arguments against this notion. First, the back- 
ground count rate is a factor of ~40 fainter than that of the 
source background: the 5-9 keV count rates in the back- 
ground and source (after background subtraction) were 
0.045 c s -1 and 1.748 c s _1 , respectively. Additionally, 
there were no differences seen in the data when screened 
against the periods of highest background rates, as dis- 
cussed in §2. Second, if it were the case that the 7.48 keV 
line is relevant in affecting the spectrum, then one might 
also expect features near 8.05 keV as well, since the 8.05 
keV line is actually stronger (Katayama et al. 2004); how- 


ever, no such features are evident in the pn spectrum here. 
Third, these instrumental lines’ fluxes are at a minimum 
within 5' of the core (Katayama et al. 2004); the back- 
ground in this case was extracted over a region located 3' 
away from the core. Fourth, we note that S05 attempted 
to detect expected instrumental features in the pn back- 
ground to verify the energy scale, but did not detect any 
instrumental lines with sufficient significance. 

Our own attempts to model these features in the source 
spectrum also argue against the notion these features were 
somehow not properly subtracted. For simplicity, we 
started with Model 8, the dual-diskline model. A third 
Gaussian component was added, with energy fixed at ei- 
ther 7.48 keV or 8.05 keV. The Gaussian width a was set to 
1 eV initially, but allowed to vary. The line normalization 
was allowed to be positive or negative. The best fits still 
yielded significant residuals near 7. 5-7. 6 keV, and failed to 
achieve as good a model fit as Model 12. We are confident 
in concluding that the feature at 7.56 keV in the observed 
frame is intrinsic to the spectrum of IC 4329a and is not 
an artifact of the instrument or the background. 

It has been reported that for pn “double event” data, 
events where photon energy is deposited on adjacent pix- 
els, the registered photon energy can be 20 eV greater than 
the energy corresponding to an otherwise-identical single- 
pixel event (e.g., Pounds & Page 2005). The pn data had 
initially been extracted using pattern 0-4 events (single 
and double events); to test for this effect, we re-analyzed 
the pn data using only single-event (pattern=0) data. For 
both filtered and unfiltered data, the pattern=0 spectrum 
was virtually identical to the pattern<4 spectrum, though 
with ~30% fewer counts, and a continuum slope that was 
steeper by ~0.05 in photon index. Importantly, however, 
the fitting results for Fe K bandpass features were virtu- 
ally identical in all cases. We conclude that the emission 
and absorption features modeled here are not artifacts of 
pattern selection. The data/model continuum residuals 
showed slight improvement above ~9 keV, but since our 
focus was quantifying the Fe K bandpass features using 
the maximum usable exposure, we continued to use the 
pattern<4 data. 

4.3.3. Estimating the significance of the detection of the 
7.68 keV absorption feature 

The standard, two-parameter F-test for the addition 
of a Gaussian to model the absorption line at 7.68 keV 
yielded a null hypothesis probability (the probability that 
the null hypothesis, that the model without an absorption 
line at 7.68 keV, is correct) of 5.2 xlO -5 (~4.1cr signifi- 
cance). However, using the F-test in this manner has a 
tendency to overestimate the detection significance, as the 
F-test does not take into account the possible range of en- 
ergies where a line might be expected to occur, nor does 
it take into account the number of bins (resolution ele- 
ments) present over that energy range. The F-test can 
yield the probability Pi (equal to one minus the null hy- 
pothesis probability) of finding a feature at a given energy 
if the line energy is known in advance (see e.g., warnings by 
Protassov et al. 2002). However, in this case, there was no 
a priori expectation of a 7.68 keV feature, and in searching 
for narrow features at arbitrary energies, one is searching 
over many resolution elements and needs to account for 
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the possibility that narrow features can occur by chance 
due to statistical noise. One might search the 4-9 keV 
bandpass for Fe K emission or absorption lines, including 
narrow features near 5-6 keV that have been generally in- 
terpreted as gravitationally redshifted Fe K emission (e.g., 
Turner et al. 2002, 2004). There are about 36 pn reso- 
lution elements over this energy range. The probability 
of detecting a feature at any energy in this range due is 
found from (Pi) N , where N is the number of resolution 
elements. For IC 4329a, the probability that the 7.68 keV 
feature is spurious thus becomes 1.9x 10 -3 . This is still 
significant at greater than 3.1a significance. It is noted 
that if the spectral bin sizes are smaller than the instru- 
ment resolution, then the number of bins where a line may 
be located has probably been over-estimated; this means 
the detection probability may have been underestimated. 

Monte Carlo simulations were performed as an addi- 
tional test of the line significance (see Porquet et al. 2005). 
These simulations tested the null hypothesis that the spec- 
trum is well-fitted by a model that does not include the 
7.68 keV absorption feature. For simplicity, Model 8 (a 
power-law plus a pexrav component, plus two disklines 
for Fe Ka and Fe f3 emission, all modified by WABS) was 
initially assumed. In order to generate 1000 fake spec- 
tra, with the photon statistics appropriate for a 93 ksec 
good time exposure, one could simply run the XSPEC com- 
mand “fakeit” on this model 1000 times. However, this 
process does not account for the uncertainty in the null 
hypothesis model itself; such uncertainty is relevant, for 
example, when one is testing for the presence of broad 
features, when broad features are present in the null hy- 
pothesis model, or when the original data set is relatively 
noisy. Instead, we initially assumed the null hypothesis 
model, ran fakeit once, with the photon statistics for a 93 
ksec exposure, and then re-fit the model to this new faked 
data set, yielding a ’’modified” null hypothesis model. We 
ran fakeit a second time, using this modified model, to 
generate a second faked spectrum. This entire process was 
repeated 1000 times to generate 1000 fake spectra (based 
on the range of the 1000 ’’modified” null hypothesis mod- 
els, as opposed to the single, original model). For each 
faked spectrum, we fit a model of the same as the null 
hypothesis model, obtaining a x 2 value (hereafter x 2 u ii)- 
We then added a narrow Gaussian component (width a 
fixed at 10 eV) to the fit. The Gaussian normalization 
was allowed to be positive or negative, to account for the 
possibility of detecting either an emission or absorption 
feature (since there had been no a priori expectation to 
expect either). We searched over the energy ranges 4.0- 
6.1 keV and 7.1-10.0 keV. The range 6. 1-7.1 keV was ex- 
cluded to avoid the possibility that the e.g., Fe K emission 
lines or 7.1 keV PEXRAV edge would bias the Monte Carlo 
results. For each faked spectrum, we stepped the Gaus- 
sian centroid energy over these energies in increments of 
0.1 keV (e.g., slightly higher than the instrument resolu- 
tion), fitting separately each time to ensure that the lowest 
value of x 2 was found. For each faked spectrum, the mini- 
mum value of x 2 was compared with x 2 u ii> yielding a Ax 2 
value. Repeating this process 1000 times yielded a distri- 
bution of 1000 Ax 2 values which could then be compared 
with the observed Ax 2 value, which for this case was - 
21.73. The result of these simulations was that a |Ax 2 | 


value of 21.73 was not exceeded even once. The values of 
| Ax 2 | corresponding to the 68, 90, 95, and 99.0% values in 
the cumulative frequency distribution were 6.8, 9.5, 11.1, 
and 14.4, respectively; the highest simulated | Ax 2 1 value 
was 18.1. The simulations thus suggest that the probabil- 
ity that the null hypothesis model was correct is at most 
1/1000, or that the 7.68 keV absorption feature is detected 
at >99.9% confidence. 

Finally, we briefly address “publication bias”: there 
are approximately a dozen XMM-Newton observations of 
Seyfert Is of similar length and quality to the current ob- 
servation. One could thus argue that the significance of de- 
tecting the 7.68 keV feature should be reduced to >98.8%, 
given that there was no a priori reason to expect an ab- 
sorption feature at such an arbitrary energy. We leave it 
to the reader to include the relevance of such publication 
bias. However, given the very high value of |Ax 2 |, the 
unlikeliness that the feature is an artifact of background 
modeling, and and the low probability that this could be 
due to photon noise, we henceforth assume that the 7.68 
keV absorption line is intrinsic to IC 4329a. 

4.3.4. The origin of the line 

As far as the origin of this feature is concerned, there 
are several candidates to consider. K-shell absorption due 
to Fe XXV or Fe XXVI is a prime candidate, given that 
other absorption lines attributed to highly-ionized Fe have 
been found in AGN. If the 7.68 keV line is due to Fe XXVI, 
then it is blueshifted by 0.72 keV, implying an origin in 
material that is blueshifted by nearly 31000 km s _1 rel- 
ative to systemic. If the line is due instead to Fe XXV, 
then the blueshift is even higher: assuming an origin in 
the resonance line (6.697 keV in the rest-frame), then a 
blueshift of 0.98 keV relative to systemic is required, sug- 
gesting an origin in matter flowing ~44000 km s -1 towards 
the observer and relative to systemic. 

To shed additional light on the origin of this feature, we 
modified the double-Gaussian model to include absorption 
due to ionized material along the fine of sight (Model 13). 
We modified Model 8 by including an XSTAR component, 
assuming solar abundances, and keeping the redshift as 
a free parameter. The best-fit model (see Table 7) was 
one with a high ionization parameter £: log(£) = 3.731° 13 
erg cm s -1 . x 2 /dof was 1809.17/1887, a slightly better fit 
compared to Model 12. The column density required was 
1.4lJ;° x 10 22 cm -2 , and the required radial velocity shift 
relative to systemic was -0.09311 qq 2 c i or ~-270001 2 q°q 
km s -1 . At this ionization level, the expected dominant 
absorption line, and the most likely origin of this feature, is 
due to Fe XXVI. The best fit model is thus consistent with 
the interpretation above that the absorption feature has 
an origin in highly- ionized material outflowing at ~0.1c 
relative to systemic. 

It is unlikely that the absorption feature is due to K/3 ab- 
sorption by highly-ionized Fe (greater than F XVII), other- 
wise one would expect strong Ka absorption by the same 
species. For example, if the 7.68 keV absorption feature is 
due to redshifted Fe K/3 XXV, with a rest frame energy 
near 7.87 keV, this implies a redshift of 0.024 relative to 
systemic. If the material which gives rise to this feature 
has an ionization parameter f of ~3.2-3.7, then we would 
expect very strong Fe Ka XXVI and Fe Ka XXV res- 


10 


Fe K Emission and Absorption in IC 4329a 


onant absorption at the same redshift, appearing at rest 
frame energies of 6.80 and 6.51 keV, respectively. Looking 
at the residuals in Models 12 or 13, there are no obvious 
residuals at these energies. Adding an inverted Gaussian 
with centroid energies fixed at either location, and width 
fixed at 35 eV (the Gaussian width of the 7.68 keV line in 
Model 12), results in no improvement to the fit, with up- 
per limits to the equivalent widths of 3 eV in either case. 
If the 7.68 keV feature is due to redshifted Fe K/3 XXVI 
(rest frame energy near 8.3 keV), the implied redshift rel- 
ative to systemic is 0.075, causing the Fe Ka XXVI and 
Fe Ka XXV lines to appear at rest frame energies of 6.48 
keV and 6.20 keV, respectively. The upper limits to ab- 
sorption at these energies are 3 eV and 8 eV, respectively, 
with no fit improvement in either case. 

Similarly, it is unlikely that the 7.68 keV feature is due 
to blueshifted K/3 absorption from lowly-ionized Fe (~7.1- 
7.2 keV rest frame). The absorbing gas would have to have 
an ionization parameter smaller than log/ ~ 2.1 in order 
for absorption near ~7. 1-7.2 keV to exist. For an absorber 
with log/=2.0 and a blueshift of 0.07c, a column density of 
1.1 xlO 22 cm -2 is required to fit the 7.68 keV absorption 
line. However, this model requires very strong spectral 
curvature below ~4 keV ; such curvature is not seen in the 
EPIC spectrum, as XSTAR modeling using low ionization 
photo-ionized absorbers yield poor fits to the pn spectrum 
below ~4 keV for column densities above roughly 10 21 
cm -2 . Additionally, for photoionized material in the ap- 
proximate range log/; — 1. 5-2.1, Fe Ka absorption is just 
as strong, or stronger than, the Fe K/3 absorption. For ma- 
terial in this range of ionization parameter, if the 7.68 keV 
line were due to neutral or lowly-ionized Fe K/3 absorp- 
tion, blueshifted by ~0.6, then we would expect neutral 
Fe Ka absorption to appear near 6.9 keV (rest frame). 
There are no obvious residuals near this energy in Models 
12 or 13. Forcing an inverted Gaussian into the model, 
with width tied to that for the 7.68 keV Gaussian line, 
yields no change in the fit, with limits to the absolute val- 
ues of the intensity and equivalent width of 2.7 x 10“ 6 ph 
cm -2 s -1 and 2.8 eV, respectively. One can alleviate the 
requirement of having Fe Ka absorption by noting that Fe 
Ka is negligible compared to Fe K/3 for very low values 
of the ionization parameter (log/' < 1.2). For these ion- 
ization parameters, the dominant species will be <XVII; 
the L-shell is full, so Ka absorption cannot be produced. 
However, as above, the required column density in this 
case would introduce very strong spectral curvature into 
the spectrum at and below ~4 keV, which is not seen. An 
origin for the 7.68 keV absorption feature in Fe K/3 thus 
seems unlikely. 

We next considered the possibility that the 7.68 keV ab- 
sorption feature could be due to very highly-blueshifted, 
moderately-ionized Fe Ka (Fe XVII-XXII or so, with 
rest-frame energies near 6.44-6.51 keV). For instance, an 
XSTAR absorber with log/=2.1±0.2, a column density Nr 
= 1.8±0.5 xlO 20 cm -2 , blueshifted by 0.16 ±0.01 c rela- 
tive to systemic provides almost as good a fit above 3 keV 
as Model 13, with x 2 /do/ = 1820.49/1887 for the best- 
fit model. Because moderately-ionized Fe has a higher 
opacity than H- or He-like Fe, a lower column density can 
account for the 7.68 keV absorption line without intoduc- 
ing strong spectral curvature around 3-4 keV. However, 


extrapolation of this model to energies down to ~0.9 keV 
reveals large data/model discrepancies. Specifically, for a 
log/=2.1 absorber with Nr ~ 2xl0 20 cm -2 , one expects 
very strong (EW = 25-30 eV) absorption features near 
rest-frame energies of ~0.91-0.94 keV due to L-shell ab- 
sorption by moderately-ionized Fe ~ XIX-XX and near 
1.86 keV, likely due to absorption in He-like Si. For an ab- 
sorber blueshifted by 0.16c, we would expect such features 
to appear at 1.06 and 2.13 keV. However, the pn data at 
those energies show no such features; upper limits to in- 
verted Gaussians fixed at those energies yield only small 
(< 1 eV) upper limits. An origin in very highly blueshifted 
Fe XVII thus seems unlikely. 

There are transitions around 7. 6-7. 8 keV associated 
with highly ionized Ni K (XXV - XXVI or so), but an 
origin in Ni K at a velocity shift that is close to or at 
systemic seems unlikely: unless Ni is many times over- 
abundant compared to Fe in IC 4329a, we would expect 
to see very strong evidence of additional absorption due 
to highly-ionized Fe Ka in the ~6.6-6.96 keV range. Even 
if the 7.68 keV feature is due to strongly redshifted He- 
or He-like Ni or strongly blueshifted low- or intermediate- 
ionized Ni, we would still expect strong evidence for ad- 
ditional Fe K absorption lines at the appropriate redshift 
or blueshift. Given that there is only one absorption line 
strongly evident in the data, an origin in Ni seems unlikely. 

A local (z— 0) origin should be discussed as well, given 
that several recent soft X-ray AGN spectra have revealed 
evidence for absorption due to the local hot interstellar 
medium (e.g., McKernan, Yaqoob & Reynolds 2004). In 
the case of IC 4329a, for this line to be local in origin, 
it would have to originate in K/3 absorption from inter- 
mediate ionized Fe (Fe ~ XXII). However, the very high 
column density required to produce such a line (>10 22 
cm -2 ) is an implausible physical description of the local 
hot interstellar medium. Similarly, the relativistic veloci- 
ties required for strongly blueshifted Fe Ka absorption are 
not associated with the local hot interstellar medium, and 
a local origin for the line is unlikely. 

Finally, Figure 7 shows a model spectrum of the effect 
of the soft X-ray absorbing components seen in S05 on a 
simple power-law continuum. These low and moderately 
ionized components are too low in column density to intro- 
duce strong spectral curvature above ~4 keV, suggesting 
that our emission and absorption modeling of the Fe K 
bandpass is robust. 

4.4. Consistency checks 

We checked for consistency between the pn and the 
MOS2 by applying the above models to the MOS2 data. 
We ignored all MOS2 data below 4.5 keV due to calibration 
uncertainties and data above 12 keV. Figure 8 shows the 
pn and MOS2 residuals to models consisting of a power-law 
plus PEXRAV, modified by neutral absorption; the shape of 
the Fe K profiles are overall similar. Applying the models 
above to the MOS2 data yielded similar and statistically 
acceptable fits. That is, most spectral parameters were 
consistent between the two instruments, though the pho- 
ton indices for the MOS2 spectra were about 0.15 steeper. 
We were able to confirm that the Fe K/3 line is detected 
in the MOS2 data at ~90% confidence, using an F-test to 
compare models 3 and 4. However, due to the lower sig- 
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nal/noise in the MOS2, it was not possible to verify e.g., 
the existence of the 7.68 keV absorption line or confirm 
that a dual-diskline model was a better fit than a double- 
Gaussian model. 

To check for consistency between the pn and the RXTE 
PCA, we generated a ’7-day’ simultaneous PCA spectrum 
using RXTE data taken from the 36 obsids which occurred 
within ±3.5 days of the midpoint of the XMM-Newton 
observation. This resulted in a PCA spectrum with a 
good exposure time of 19.6 ksec, enough to confirm that 
an Fe Ka line is present, but not enough to detect other 
Fe K bandpass spectral features. The photon index for 
Model 3, applied to the ’7-day’ spectrum, was 1.89±0.4, 
steeper than that of the model when applied to the pn 
data. The 2-10 keV flux from the ’7-day’, medium-term 
and long-term modeled spectra are ll.OxlO -11 erg cm -2 
s -1 , 11.2xl0 -11 erg cm -2 s -1 , and 12.0xl0 -11 erg cm -2 
s -1 , respectively. These values are about 30% than the 2- 
10 keV flux estimated from model fits to the XMM-Newton 
pn spectrum, 8.9xl0 -11 erg cm -2 s -1 . This discrepancy 
with the PCA normalization has been noted previously; 
see e.g., Yaqoob et al. (2003) 3 . 

5. TIME-RESOLVED SPECTRAL FITTING 

Studying the time-resolved variability of the Fe K line 
and continuum components provides a complementary 
analysis to single-epoch spectral fitting. The simplest 
models predict that if the reprocessing time scale is negligi- 
ble, then variations in the Fe K line flux should track con- 
tinuum variations, modified only by a time delay equal to 
the light travel time between the X-ray continuum source 
and the Fe K line origin. However, recent studies of both 
broad and narrow Fe K lines have not been able to sup- 
port this picture (e.g., Iwasawaet al. 1996, 1999; Reynolds 
et al. 2000; Vaughan & Edelson 2001, Markowitz, Edel- 
son & Vaughan 2003). Most of these studies found that 
Fe K lines tend to vary much less than the continuum, 
with no evidence for correlated variability, on any time 
scale. Markowitz, Edelson & Vaughan (2003) addition- 
ally showed that Fe K lines, like the continuum, tend to 
exhibit stronger variability amplitudes towards relatively 
longer time scales. 

Variability of narrow absorption features has also been 
reported, e.g., Reeves et al. (2004) found evidence for the 
narrow Fe XXV absorption line in the XMM-Newton spec- 
trum of NGC 3783 to increase in equivalent width over ~a 
day, suggesting an origin close to the central X-ray illumi- 
nating source. We investigated time-resolved spectroscopy 
IC 4329a to search for changes in the profile, fluxes, or 
peak energy of the Fe K emission or 7.68 keV absorption 
features in IC 4329a, and explore if any such variations are 
linked to continuum flux variations. 

We started with short-term data and split the pn data 
into two time halves, each covering a duration of 66.3 ksec 
and with a good time exposure of 46.4 ksec. We fit each 
spectrum over 2.5-12 keV with Model 12 (dual-diskline 
plus narrow Gaussian in absorption). The best-fit spec- 
tral parameters, listed in Table 8, are all consistent with 
each other; differences in best-fit values can likely be at- 
tributed to less-than-optimal photon statistics. There is no 
evidence for temporal evolution in any of the absorption 


and emission profile parameters. The data/model residu- 
als when each spectrum is fitted with a simple absorbed 
power-law are plotted in Figure 9 and show no obvious 
changes in the profile. However, given that the hard X- 
ray continuum light curve (shown in S05) displays mini- 
mal variability on these time scales, this is not surprising. 
Additionally, given the minimal continuum variability, we 
did not investigate spectral variability as a function of flux, 
i.e., high/low-flux states. 

Next, we performed high time-resolution spectroscopy 
on all time scales to yield light curves of T and the Fe K 
flux. Time bins were chosen with the goals of achieving 
adequate signal-to-noise in the Fe K flux light curve. The 
medium-term data were divided into 18 bins, each cov- 
ering a duration of 1.972 days, or once every 30 satellite 
orbits; each bin contained ~10 separate RATE visits. The 
average exposure time per bin was 5.6 ksec. The long-term 
data were divided into 18 bins of approximately 44 days 
each in duration, with no bin overlapping the two 56-day 
monitoring gaps: data from the periods 2003 Apr 08 - 2003 
Sep 29, 2003 Nov 24 - 2004 Sep 30, and 2004 Nov 25 - 2005 
Oct 02 were divided into 4, 7, and 7 bins, respectively. The 
average good time per bin was 8.1 ksec, excluding bins 3 
(100.9 ksec) and 4 (68.7 ksec), which overlapped with the 
medium-term sampling. 

Models were applied to time-resolved data based on the 
time-averaged fits. The short-term data were divided into 
8 time bins, with each bin covering a duration of 16.6 ksec 
a with a good time exposure of 11.6 ksec. Spectral fit- 
ting was done over 2.5-12 keV. Model 5 (see Table 2), 
which featured Gaussians for the Fe Ka and K/3 emission 
lines, was used for the short-term spectra. Using the dual- 
diskline model offered no improvement to the fits; also, 
the 7.68 keV absorption line was not detected in any time- 
resolved spectrum. For both the long- and medium-term 
data, Model 3 (see Table 1), which featured a single Gaus- 
sian to model the Fe Ka emission, was used. More de- 
tailed fits were not practical, e.g., the Fe K/3 line was not 
detected. Spectral fitting was done over 3.5-24 keV. The 
HEXTE statistical errors were too large to study variabil- 
ity, so HEXTE was excluded from the time-resolved anal- 
ysis. The Fe K a line detected in nearly every time bin at 
99.8% confidence or greater in an F-test. When fitting the 
time-average models, it was not practical to thaw param- 
eters such as the Fe K line profile shape or peak energy, 
given the PCA resolution and sensitivity. The neutral ab- 
sorbing column density was kept as a free parameter. The 
F-test was used to determine which free parameters to 
include in the fits. The time-resolved spectral fits were 
repeated with each of photon index T, reflection fraction 
R, and Fe Ka line normalization Ipe frozen at their time- 
averaged values, and the values of total y 2 were compared 
via an F-test to determine which parameters could justifi- 
ably allowed to be free in the fits. The results are shown 
in Table 9. We found no formal requirement to thaw R on 
any time scale: the total values of x 2 decreased only by 
0.06 or 0.01 (long and medium, respectively) when R was 
thawed from the time-averaged value. All analysis here- 
after assumes R is kept fixed at the time-averaged value 
for each time scale. As shown in Table 9, it was significant 
at >99.99 confidence to thaw T only on the medium and 


3 This issue is also discussed in an ASCA Guest Observer Facility Calibration memo at http://heasarc.gsfc.nasa.gov/docs/asca/calibration/3c273.results 
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long time scales, suggesting statistically significant varia- 
tions in T on these time scales. It was significant at >99% 
confidence to thaw Iy e only on the long time scale. Errors 
on the Fe K line flux and T light curves were determined 
using the point-to-point variance, as discussed in detail 
by Vaughan & Edelson (2001) and Markowitz, Edelson & 
Vaughan (2003). The 2-10 keV flux F 2 _io was measured 
from each model fit. Errors on F 2 _io within a time bin 
were derived from the 2-10 keV continuum light curves, 
using the distribution of points within that time bin. 

Mean spectral fit values and errors are listed in Table 
10. Figure 10 shows the light curves for F 2 _io, F, and 
/ Fe for all three time scales We note that the long-term 
averages to the spectral fit parameters shown in Figure 10 
and listed in Table 10 do precisely match the spectral-fit 
parameters to the long-term time average-spectrum listed 
in Table 1. (see §4.4 for notes comparing the PC A and pn 
2-10 keV relative flux normalizations). This is due to the 
fact that the long-term time-average PCA spectrum has 
a substantial contribution from the medium-term time- 
average spectrum (the medium-term time-average spec- 
trum comprises 36% of the long-term spectrum’s total ex- 
posure time). Figure 11 shows correlation diagrams for 
r and Jpe plotted against F^-io- Table 10 also lists the 
Pearson correlation coefficients and null hypothesis proba- 
bilities. Finally, the fractional variability amplitudes, F var , 
as defined in Vaughan et al. (2003), are listed in Table 11 
for the light curves. 

From Figures 10 and 11 and Tables 10 and 11, one can 
see that there is not much variability in F 2-10 or F on 
short time scales. However, F tends to correlate well with 
F 2 _io on all three time scales, despite the somewhat nar- 
row range spanned by F 2 _io. However, there is no strong 
evidence for variability in Ip e on any time scale; null hy- 
pothesis probabilities for a constant line flux cannot be 
rejected at greater than 60% confidence on any time scale. 
The lack of variability of the Fe line flux on short time 
scales has already been reported by S05. The results on 
long time scales are similar to those of Weaver, Gelbord & 
Yaqoob (2001), whose did not find any evidence for strong 
Fe line flux variability in IC 4329a based on five separate 
ASCA observations between 1993 and 1997. 

6. DISCUSSION 

6.1. The Fe K emission profile: moderately-relativistic 

disklines 

The XMM-Newton pn spectrum of IC 4329a reveals an 
Fe K emission profile with two peaks, identified as Fe Ka at 
6.4 keV and Fe K/3 at 7.0 peak. The 6.4 keV peak is much 
stronger in terms of both EW and height above the power- 
law continuum than the 7.0 keV peak, ruling out models 
wherein both peaks are the Doppler horns of a single rel- 
ativistic diskline model (e.g., McKernan & Yaqoob 2004). 
The red peak in particular is not symmetric: a moderate 
red wing is present, and models incorporating two symmet- 
ric peaks do not provide as good a fit to the data as mod- 
els incorporating at least moderately-relativistic diskline 
emission from Fe Ka and Fe K/3. IC 4329a’s Fe K emis- 
sion profile is thus distinct from more symmetric and more 
narrow profiles which lack a red wing, such as that seen 
in NGC 5548 (Yaqoob et al. 2001, Yaqoob & Padmanab- 
han 2004). Such narrow, symmetric profiles are generally 


suspected of having an origin in matter far from the black 
hole, such as in the outer accretion disk, broad line region, 
or the molecular torus. In IC 4329a, however, the presence 
of a moderate red wing leads us to conclude that the emis- 
sion profile is instead modeled equally well with either of 
the following two scenarios: 1) each peak is modeled by a 
low-inclination, moderately-broadened diskline, plus little 
or no contribution from a narrow emission line, or 2) each 
peak is modeled as the sum of a narrow line, to account 
for the emission closest to the rest energy of the line, plus 
a high-inclination diskline. In the latter case, neither the 
narrow nor broad component dominates the emission. 

In some Seyferts, identification of broad Fe lines is com- 
plicated by model degeneracies which arise when 3-6 keV 
spectral curvature due to warm absorption is taken into 
account (e.g., Reeves et al. 2004, Turner et al. 2005). In 
IC‘4329a, however, only minimal spectral curvature above 
3 keV, and no spectral curvature above 5 keV, is expected 
given the results of soft X-ray absorption modeling (see 
Figure 7), and so the detection of a moderate red wing 
is robust. The best fitting model includes emission ex- 
tending in to a radius of 13jjJj 8 R g , similar to that found 
by Done et al. (2000) using ASCA. Done et al. (2000) 
note that IC 4329a’s profile does not seem as broadened 
and reddened as that of MCG-6-30-15, where the diskline 
emission extends down to a ~2-3 R g (Wilms et al. 2001; 
Young et al. 2005). The red wing emission in IC 4329a 
extends down to energies of ~5.4 keV (Done et al. 2000) 
or ~5.7 keV (this work), as opposed to energies below 4-5 
keV in the case of MCG-6-30-15. Although our analy- 
sis cannot formally rule out an inner emission radius of 
<2f? g , the Fe K emission profile in IC 4329a thus seems 
to fall into a category between those consisting of narrow 
lines and those consisting of very strongly relativistically 
broadened lines. 

Assuming that the bulk of the Fe K emission originates 
in the inner regions of an optically-thick, geometrically- 
thin accretion disk, as in case 1) above, these findings 
are consistent with the notion that the accretion disk in 
IC 4329a could be truncated below ~10-20 R g , or could 
transition from a geometrically- thin, optically-thick disk 
at large radii to a radiatively-inefficient flow, such as an 
advection-dominated accretion flow (ADAF; Narayan & 
Yi 1995), in the innermost regions of the accretion disk, 
e.g., Dove et al. (1997) and Esin, McClintock & Narayan 
(1997). 

Additionally, for case 1), the best-fit rest energy for 
the Fe Ka diskline was 6.44±0.01 keV, more consistent 
with an origin in moderately-ionized Fe than in neutral 
Fe; the likely dominant ionization stage is ~XVII. This 
phenomenon could be related to IC 4329a’s high X-ray 
luminosity. Using the XMM-Newton 2-10 keV flux of 
9x 10 -11 erg cm -2 s -1 , we derive an 2-10 keV luminos- 
ity L 2 -io of 5.1 xlO 43 erg s _1 (assuming H a = 70 km s -1 
Mpc -1 and A o =0.73). Similarly, Yaqoob & Padmanabhan 
(2004) measured Fe line energies that seemed more con- 
sistent with moderately-ionized Fe than neutral Fe in two 
other high-luminosity Seyferts, 3C 120 and Mkn 509. 

The lack of variability in the Fe K emission line is con- 
sistent with an origin in distant matter, as suggested e.g., 
by S05. However, the line profile suggests a non-negligible 
contribution from material origin close to the black hole. 
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If in fact the bulk of the Fe K emission originates in the 
inner accretion disk, as in case 1), then a likely scenario is 
as follows: assuming that variations in the line track those 
in the X-ray continuum, modified only by a light travel 
time delay, the line will not vary strongly, unless the con- 
tinuum does. However, the continuum in IC 4329a does 
not display strong variability. Even on time scales of ~2 
years, the peak-to-trough continuum variations do not ex- 
ceed ~50%, which simply may not be sufficient to trigger 
very strong trends in Fe K flux. Similarly, with time bins 
chosen to maximize the variability-to-noise in the in the 
Fe flux light curves, this analysis may simply not be suf- 
ficiently sensitive to track small-amplitude trends (<50% 
variations) in line flux. 

If there are approximately equal broad and narrow com- 
ponent contributions to the profile, as in case 2), then, 
for time scales longer than the light-crossing times be- 
tween the continuum source and the sites of broad and 
narrow line production, it could again be the case that the 
lack of variability in either component is a result of the 
lack of strong continuum variability. However, the PCA is 
not sensitive enough to distinguish between the broad and 
narrow components. Therefore, if the narrow component, 
which comprises ~60% of the total (narrow plus broad) is 
constant, then the PCA time-resolved spectroscopy may 
not be sensitive to small changes in the broad component 
only on time scales longer than a few days. In any case, 
the lack of strongly-correlated variations between the con- 
tinuum and line prevents us from making any firms con- 
clusions about using variability to determine the location 
of the line-emitting material. 

6.2. The 7.68 keV absorption feature: the highly-ionized 
outflow in IC 4329a 

High-resolution spectroscopy of the Fe K band of sev- 
eral AGN has yielded evidence for highly-ionized absorbing 
gas, e.g., NGC 3783 (Reeves et al. 2004) and MCG-6-30-15 
(Young et al. 2005). In addition, there is evidence in PG 
and BAL quasars for Fe K lines and/or edges which suggest 
absorbing gas which is both highly-ionized and outflowing 
at relativistic velocities. For instance, in PG 1211+143 
(Pounds et al. 2003a), PG 0844+349 (Pounds et al. 2003b), 
PDS 456 (Reeves et al.2003), and the BAL quasars APM 
08279+5255 (Chartas et al. 2002; Hasinger, Schartel & Ko- 
mossa 2002), and PG 1115+080 (Chartas et al. 2003), the 
ionization parameters were typically log£ ~ 2.5 - 3.7, the 
outflow velocities spanned 0.08-0. 34c, and the absorber 
column densities were typically ~l-5xl0 23 cm -2 , (al- 
though see McKernan, Yaqoob Sz Reynolds 2004 regarding 
caveats for PDS 456 and PG 1211+143). 

The highly-ionized, high-velocity absorbing material in 
IC 4329a has a derived ionization parameter and outflow 
velocity similar to what is found in these objects, though 
we note that the derived column density for the IC 4329a’s 
outflow component is about a factor of ten lower. Nonethe- 
less, we hereby add IC 4329a to this rather small but grow- 
ing list, noting that IC 4329a is thus the lowest-redshift 
AGN known to exhibit this type of outflow. 

To estimate the distance r between the central black 
hole and the outflowing gas, we can use £ = L 2-200 / ( n 
r 2 ), where n is the number density. L 2-200 is the lumi- 
nosity corresponding to the 2-200 keV illuminating contin- 


uum, which we use instead of the 1-1000 Rydberg ionizing 
continuum luminosity because most of the photons that 
ionize this particular component are expected to be hard, 
not soft, X-rays. We estimate the maximum possible dis- 
tance to the material by assuming that the the thickness 
A r must be less than less than the distance r. Noting that 
the column density Nh = nAr yields the upper limit r < 
L 2 - 20 o/(N h £). We estimate the 2-200 keV flux from our 
models to be 4.9xlO -10 erg cm -2 s -1 , which corresponds 
to L 2-200 = 2.8xl0 44 erg s _1 . This yields r < 3.8xl0 18 
cm, or 1.2 pc. 

In the forementioned PG and BAL quasars, the derived 
mass outflow rates are usually at least 1 M© yr _1 and 
on the order of the accretion inflow rate, suggesting that 
these outflows account for a substantial portion of the to- 
tal energy budget of the AGN. Under the assumption that 
the gas is in equilibrium and that the outflow velocity is 
a constant, the mass outflow rate M out for the IC 4329a 
component can be derived via conservation of mass: M out 
= flnr 2 vm p , where v is the outflow velocity, m p is the 
proton mass, and fl is the covering fraction, which we will 
assume is 47r/10 sr. Substituting nr 2 = Z, 2 - 2 oo/£, we And 
an outflow rate of ~ 3xl0 27 gm s -1 = 5 M© yr _1 . We 
note that the actual outflow rate should be lower if there 
is an extreme degree of collimation, or higher if the out- 
flowing component is not directed along the line of sight. 

We can compare this to the inflow accretion rate M acc 
using Lboi = f?M acc c 2 , where t] is the accretion efficiency 
parameter, typically 0.1. Using the relation of Padovani & 
Rafanelli (1988), Lboi = 59 vL v at i/=2 keV, and a photon 
index T of 1.75, the bolometric luminosity is thus esti- 
mated to be Lboi = 321,2- 10 = 1.6xl0 45 erg s -1 . We find 
M acc ~ 0.3 M© yr -1 . The kinetic power associated with 
the outflow component, estimated as M ou t+ 2 , is 3xl0 45 erg 
s -1 , a factor of ~60 higher than the bolometric luminos- 
ity. Similar to the situation in the aforementioned quasars, 
the outflow rate of this particular absorbing component 
may therefore represent a substantial fraction of the ac- 
cretion inflow used to power the bolometric luminosity in 
IC 4329a. Similar calculations on the low-velocity, soft 
X-ray absorbing components (S05), using Li on , estimated 
to be 1.4 xlO 44 erg s _1 , suggest that they, too, are associ- 
ated with outflow rates that are of the same magnitude or 
higher than the inflow accretion rate. 

The high-velocities of these outflows suggests that the 
outflow is likely to be associated with the accretion disk at 
a small radius. An explanation commonly invoked for the 
high-velocity, high-ionization absorbers in that the flow 
originates in the innermost part of a radiatively-driven ac- 
cretion disk wind (e.g., Proga, Stone & Kallman 2002; 
King & Pounds 2003). In this scenario, the winds must 
be launched from very close in to the black hole, yet far 
enough away from the black hole that the outflow velocity 
exceeds the local escape velocity. For a Keplerian disk, the 
radius at which the escape velocity is 0.1c will be equal to 
(v esc /c) 2 Rsch = lOORsch! a radiatively-driven wind could 
thus be applicable to IC 4329a if (neglecting acceleration) 
the wind is launched from a larger radius. 

More specifically, King & Pounds (2003; see also Reeves 
et al. 2003) demonstrated that objects accreting near 
the Eddington rate are likely to exhibit radiatively-driven 
winds with substantial column densities, and the outflow 
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rates of these winds may be comparable to the accretion in- 
flow rate. If the flow becomes optically-thick at small radii 
(~10-100 R g ), extreme velocities and high column densi- 
ties are likely to the associated with the flow. This sce- 
nario may be applicable to the high-velocity absorber in IC 
4329a. However, we do not know a priori if the accretion 
rate of IC 4329a is indeed near Eddington. Uncertainty re- 
garding IC 4329a’s black hole mass Mbh prevents precise 
knowledge about the accretion rate. The reverberation- 
mapped estimate for .Mbh (Peterson et al. 2004) is for- 
mally an upper limit only, 9.9+g 7 9 9 M e , implying an accre- 
tion rate relative to Eddington of Z/boi/^Edd of 1.29lo n 8 f 3 . 
Other methods yield higher black hole mass estimates. For 
example, using the photo-ionization method (estimating 
velocity dispersions based on H/3 line widths) to estimate 
the distance from the central illuminating source to the 
broad-line region, Wandel, Peterson & Malkan (1999) es- 
timate a black hole mass of 2.2xl0 7 M©. Along another 
track, Nikolajuk, Papadakis k, Czerny (2004) suggest a 
prescription to estimate Mbh based on short-term X-ray 
variability amplitude measurement, specifically the nor- 
malized excess variance, ctnxs- They estimate Mbh = 
1.24x10® M© based on short-term RXTE light curves. 
We rebinned the short-term, pn 2.5-12 keV light curve 
with a sampling time of 6.8 ksec, to yield a light curve with 
enough data points (20) to get an accurate measurement of 
o’nxs- We found <tnxs= 0.00136 ± 0.00010 (errors derived 
using Vaughan et al. 2003); the Nikolajuk et al. (2004) 
relation yields a Mbh estimate of 8.64F0.60 x 10' M©, 
similar to the Nikolajuk et al. (2004) estimate 4 * . This value 
of the black hole mass places the estimated accretion rate 
Tbol/^Edd at 15±1%. A more accurate black hole mass 
determination, e.g., via more accurate reverberation map- 
ping, is needed to resolve this issue. We can only conclude 
that we currently cannot rule out an accretion rate close 
to the Eddington limit, and so a radiatively-driven disk 
wind cannot be ruled out. 

Another problem with the disk-wind scenario is that 
the 7.68 keV absorption feature is narrow, and consistent 
with a single component only. A broad range of velocities 
seems more physically plausible if we are seeing a con- 
tinuous stream of gas being accelerated from rest to 0.1c 
(relative to the systemic velocity). However, there are no 
such indications in the Fe K bandpass of material at lower 
velocities. These facts could suggest that we are witness- 
ing a discrete, transient, outflowing ’blob’ of material as 
opposed, to a continuous flow. In this case, the average 
mass outflow rate will be lower than the estimate given 
above, depending on the duty cycle of ejection. 

In addition to radiatively-driven outflows, another 
mechanism whereby radio-quiet AGN can have material 
outflowing at relativistic velocities is described by the so- 
called ” aborted-jet model” of Ghisellini, Haardt & Matt 
(2004). In this model, radio-quiet accreting black holes 
launch jets only intermittently, yielding discrete, fast- 
moving blobs that travel along the black hole rotation axis, 
as opposed to a jet that is continuously on. Extraction of 


black hole rotational energy, in addition to accretion en- 
ergy, provides the energy source. The jet is launched at 
a velocity less the escape velocity, so blobs reach a maxi- 
mum radius from the black hole before falling backwards 
and colliding with later-produced outflowing blobs. In the 
case of IC 4329a, we could thus be witnessing a singular, 
transient blob at a particular point in its outward journey, 
although its location must be < 100 Rs c h for this scenario 
to apply. 

However, even if the gas is in the form of a discrete 
blob, this particular model may not be relevant. We can- 
not ascertain whether or not the material lies on the jet 
axis, and, due to uncertainty regarding the location of the 
material, it is difficult to ascertain whether the escape ve- 
locity exceeds the material’s outflow velocity. Finally, the 
Fe K bandpass spectrum of IC 4329a yields no evidence for 
other discrete, absorbing blobs at any other velocity. Sim- 
ilarly, there are no indications from the soft X-rays (e.g., 
S05) for absorption at a similar relativistic outflow veloc- 
ity, so the absorbing material would have to be in the form 
of a single blob consisting of a high-ionization component 
only, with no lowly- or moderately-ionized contribution, 
for such a scenario to be applicable. 

Finally, we note that this high-velocity hard X-ray ab- 
sorber is unlikely to be directly physically connected to 
the soft X-ray absorbers seen by S05. The soft X-ray 
absorbers are likely too low in outflow velocity, column 
density and ionization state to be part of the same out- 
flow phenomenon as the hard X-ray absorber, suggest- 
ing physically distinct phenomena. There is no evidence 
that the soft X-ray absorbing material “feeds” or is “fed 
by” the hard X-ray absorber, and no obvious accelera- 
tion/deceleration mechanism. 

7. CONCLUSIONS 

We have presented a re-analysis of the XMM-Newton 
pn long-look spectrum of IC 4329a. The Fe K bandpass 
is dominated by two peaks, identified as Fe Ka and K/3 
emission. The 6.4 keV peak is much higher in equivalent 
width and height above the continuum compared to the 
7.0 keV peak, ruling out models whereby both peaks are 
the Doppler horns of a single relativistic diskline. Mod- 
els featuring dual diskline components provide a better 
fit than a simple dual-Gaussian model, and demonstrate 
the necessity of including diskline components in emission 
models. However, there is a degeneracy present, in that 
the emission profile can be described equally well by ei- 
ther a model in which low-inclination dual-diskline profiles 
dominate the peaks, with little or no contribution from a 
narrow component, or a model in which approximately 
equal contributions from a high-inclination diskline and a 
narrow component describe each peak. The resolution of 
an X-ray calorimeter is needed to break the degeneracy 
between these two cases and deconvolve the narrow and 
broad components. There is no strong evidence for emis- 
sion or absorption at the systemic velocity by Fe XXV or 
Fe XXVI; there is no strong evidence for a Compton shoul- 


4 We note that the Nikolajuk et al. (2004) method assumes that the power spectral density function (PSD) contains a ’break’ (change in 

power-law slope from —2 to —1) at temporal frequencies corresponding to time scales larger than the duration over which one is measuring 

&NXS (see Nikolajuk et al. 2004 for derivation, as well as additional assumptions and caveats). PDS measurement for IC 4329a is still in 
progress, pending accumulation of long-term monitoring data, (Markowitz et al., in prep.), but we will assume that the PSD break for IC 4329a 
corresponds a time scale longer than 1.5 days (as expected for an AGN with a black hole mass larger than 5xl0 6 Mg; see Markowitz et al. 
2003). 
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der. In addition, we have used RXTE monitoring data to 
extract a PCA + HEXTE spectrum to derive tight con- 
straints on the strength of the reflection component. 

We have performed time-resolved spectral fitting on the 
XMM-Newton long-look and the RXTE monitoring data 
to probe the variability of the X-ray continuum, photon 
index and Fe K line flux on timescales spanning <a day 
to 2 years. The photon index and 2-10 keV flux are well- 
correlated. However, there is no strong evidence for vari- 
ability in the Fe K line on any time scale probed, likely 
due to the minimal level of continuum variability. 

We find strong evidence for a narrow absorption line 
at 7.68 keV. Monte Carlo simulations are used to verify 
the high level of detection significance. The most likely 
origin for the line is Fe XXVI, in highly-ionized mate- 
rial blueshifted by ~0.1c relative to systemic. This phe- 
nomenon is similar to what is seen in other high lumi- 
nosity AGN, including PG and BAL quasars, making IC 
4329a the lowest redshift AGN known to exhibit this phe- 
nomenon. The mass outflow rate of this component is 
estimated to be larger than the inflow accretion rate, sug- 
gesting that the outflow represents at least a substantial 
portion of the overall energy budget. The outflow may be 
associated with a radiatively-driven disk wind. However, 


because of the narrow velocity range seen, it is also possi- 
ble that the absorbing material is in the form of a discrete 
’blob’ of emission. It is important to find as many of these 
absorbing features as possible, to guage the frequency of 
occurrence in Seyfert AGN. 
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Table 1 

Joint PCA/HEXTE fitting results for Model 3: power-law + pexrav + Gaussian 


Parameter 

Long-term 

Med-term 

X 2 /dof 

208.14/123 

82.7/101 

2-10 keV flux (erg cm 2 s _1 ) 

1.20X10- 10 

1.12x 10 -10 

25-100 keV flux (erg cm -2 s -1 ) 

1.78x lO -10 

1.74X10- 10 

r 

l.oyfL__ 0 01 5 

l-782±g;gj? 

R 

0.51±0.04 

0.35i°;°® 

E (keV) 

R qq"I"0.05 

*>•40—0 .02 

6.33 ± 0.03 

a (eV) 

228 ± 50 

341 ± 54 

/ (10 -5 ph cm -2 s -1 ) 

18.2+5'.i 

25.3A1.9 

EW (eV) 

128±l° 

197 ± 15 


Note. — Model parameters for the best-fitting Model 3 (power-law + PEXRAV + Gaussian) to the PCA/HEXTE data. Errors 
are 90% confidence for one interesting parameter. E, a, K , and EW are the centroid energy, width, intensity and equivalent 
width for the Fe Kq Gaussian emission component. Fe K/3 line emission was not considered here. The high-energy cutoff in 
the PEXRAV component was fixed at 270 keV (e.g., Perola et al. 1999); the inclination was kept fixed at 30 degrees. The 25-100 
keV flux is based on average of HEXTE cluster A and B model fluxes. 


Table 2 


Initial EPIC-pn model fits for Fe K emission lines 




Model 1 

Model 2 

Model 3 

Model 4 

Model 5 

X* /do f 


2822.59/1896 

1945.40/1893 

1865.73/1893 

1844.88/1890 

1852.78/1893 

r 


l.b5D_o.oil 

1 <^*71 + 0.012 
i -0.015 

1.737± 0.012 

1.743± 0.013 

1.742± 0.013 

Red peak 

E (keV) 

6.40 ± 0.01 

6.40 ± 0.01 

6.39 ± 0.01 

6.40 ± 0.01 


<71 


106 ± 14 

91 ± 15 

91 ± 13 

93+l| 


Ik<* 


8.8 ± 0.7 

7.a±g;S 

8.0^;? 

8.0±0.6 


EW Ka (eV) 


8 6±6 

75-tt 


77±6 

Blue peak 

E (keV) 




6-93+QQg 

7.06 (fixed) 


<72 




121 til 

(=«7l) 


I K/3 




i-3t8:S 

(= 0.11 * I Ka ) 


EWyip (eV) 




15ls 

(= 0.11 * EW Ka ) 


Note. — Line intensities are in units of 10 5 ph cm 2 s 1 . A PEXRAV component was included in all models except 1 and 
2, with reflection strength R fixed at 0.51. 2-10 keV flux for Models 4-5 was 8.92xl0 -11 erg cm -2 s -1 . 


Table 3 


EPIC-pn model fits: Attempts to model Compton shoulder emission 




Model 6A 

Model 6B 

X 2 / dof 


1839.03/1889 

1833.68/1888 

r 


1.743 ± 0.012 

1.744 ± 0.012 

Red peak 

E (keV) 

6.41 ± 0.01 

6.40 ± 0.01 


71 

71+17 

' 1 — 14 

7^+24 

‘°-19 


bfa (10 -5 ph cm -2 s- 1 ) 

6.7A0.9 

7 4+1-1 


EW Ka (eV) 

63±8 

70l“ 

Blue peak 

E (keV) 

7.06 (fixed) 

7.06 (fixed) 



(=<7l) 

(=<7l) 


1k(3 (10 -5 ph cm -2 s- 1 ) 

(= 0.11 * I Ka ) 

(= 0.11 * I Ka ) 


EW k/3 (eV) 

(= 0.11 * EW Ka ) 

(= 0.11 * EW Ka ) 

Compton shoulder 

E (keV) 

6.24 (fixed) 

6 13+ 0 - 07 
l" 3 — 0.25 


02 

(=<7l) 

(=71 ) 


I (10 -5 ph cm' 2 s -1 ) 

1.4A0.8 

1.1A0.5 


EW (eV) 

11±8 

9±5 


Note. — A pexrav component was included, with reflection strength R fixed at 0.51. A comparison of these models with the 
dual-diskline model suggests that the excess residual emission on the red side of the Ka core is emission due to the red wing of 
a moderately broadened diskline profile, and not a Compton shoulder. 
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Table 4 

EPIC-pn model fits using a single diskline component (Model 7) 


X 2 / do f 

1922.56/1893 

r 

1.739 ± 0.013 

E (keV) 

6.40 ± 0.01 

I (10 -5 ph cm -2 s _1 ) 


EW 


P 

6 (unconstrained) 

Rin 

300^270 

Inch 

30° (fixed) 


Note. - — For the Laor diskline component, the outer radius was kept fixed at 400 R s . Results are presented here for an 
inclination angle of 30°; angles of 50° or 70° yielded significantly worse fits. A PEXRAV component with R=0.51 was included. 


Table 5 

EPIC-pn model fits using dual diskline components (Model 8) 


X 2 / do f 


1835.23/1890 

r 


1.741 ± 0.013 

Red peak 

E (keV) 

6.44 ± 0.01 


I (10~ 5 ph cm -2 s _1 ) 

8 9 +1 - 7 
-0.8 


EW 

86+^ 7 

Blue peak 

E (keV) 

6.98 ± 0.09 

P 


2.1 ± 0.3 

R\n 


12+13 

Inch 


l+llo 

1 -1 


Note. — For each Laor diskline component, the outer radius was kept fixed at 400 R e . The normalization of the blue peak 
diskline was kept fixed at 0.11 times that of the red peak. A PEXRAV component with f£=0.51 was included. 


Table 6 

EPIC-pn model fits using dual-disklines and dual-gaussians (Models 9A and 9B) 





Model 9A 

Model 9B 

X 2 /dof 



1833.82 / 1889 

1835.44 / 1888 

V 



1.741 ± 0.013 

1.743 ± 0.013 

Gaussians 

Red peak 

E (keV) 

6.40 (fixed) 

6.39 ± 0.01 



I Ka, narrow 

2-5 tit 

e-iisA 



^^Kq, narrow (®V) 

oo+30 

^-22 

sei S 


Blue peak 

E (keV) 

7.06 (fixed) 

7.06 (fixed) 


(7 


50 (fixed) 

66±g 

Laor Disklines 

Red peak 

E (keV) 

6 44 + 0 - 19 
D '" — 0.15 

6-301®;" 



^Ka, broad 

broad (eV) 

6.3±3.0 

4.0±1.8 



59±28 

37±17 


Blue peak 

E (keV) 

6-96tS'.4? 

7.06 (fixed) 


P 


2-61J1 



Rin 


iq+45 
1 — 18.8 

i 9+80 


i 


11+ 1 9 
■ Li ~ll 

4 c + 19 0 
4 °-41 


Note. — Best-fit spectral parameters for a model in which the red peak (Fe Ka) and blue peak (Fe K/3) are each modeled 
by the sum of a narrow Gaussian and a broad Laor diskline. The widths cr of the two Gaussians were tied. The intensity of 
the blue peak Gaussian was kept fixed at 17/150 that of the red peak Gaussian. The inner radii J?j n , emissivity index (3 and 
inclination i of the two diskline components were kept tied. The outer radii of both disklines were kept fixed at 400 R s . Line 
intensities are in units of 10~ 5 ph cm -2 s -1 . The intensity of the blue peak diskline was kept fixed at 17/150 that of the red 
peak diskline. To get a reasonable fit, the centroid energy of the K/3 Gaussian and the rest energy of the K/3 diskline were both 
kept fixed at 7.06 keV. A PEXRAV component with R=0.51 was included in the model. The 2-10 keV flux for the best-fit model 
was 8.9xl0 -11 erg cm -2 s -1 . 
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Table 7 


EPIC-pn model fits for 7.68 keV absorption feature 




Model 10 

Model 11 

Model 12 

Model 13 

x 2 /do J 


1829.70/1889 

1819.20/1898 

1813.50/1887 

1809.17/1887 

r 


1 746"*” 0 ' 013 
I. <'40_q oi 2 

1 71 e:+0.014 
1 . 1 io_o 006 

1.734± 0.012 

1 '7OQ+0.014 
1 . i ^°_o.009 

Red peak 

E (keV) 

6.44±0.01 

6.44±0.01 

6.44±0.01 

6.44±0.01 

/ ( 10 -5 ph cm -2 s -1 ) 

o 7+0.9 
- 0.8 

8.0±0.7 

o e+0.9 
«. 0_ 0 7 

O q + 0.9 

O-*J_ 0.6 


EW (eV) 

83^® 

76±6 

8 ii® 

79 ±1 

Blue peak 

E (keV) 

6-98i°;?{ 

O QQ+ 0 - 1 1 

b.9o_o 10 

6.98±0.09 

6.98±0.09 

0 


2 -I +0.4 
0.3 

2 i + 0.2 
z - i - 0.6 

2 l +0 - 4 
z * a -0.3 

2 l +0 - 4 
z * 1 -0.3 



1 o + 13 
11-8 

1 7+26 
A —15.8 

15+ 17 

10 -13.8 

i0 -14.8 

Inch 


q+9o 

**-3 

r»+ 10 o 

U -0 

r »+ 12 o 
U -0 

rv+ 10 o 

U -0 

PEXRIV 

Edge 

log? (erg cm s _1 ) 
E (keV) 

<0.04 

7 9q 0.25 
'• zy -0.17 



optical depth r 


0.034 ± 0.013 

7 - 68 iao3 


Gaussian 

E (keV) 





a 



35±|1 



| I | ( 10 -5 ph cm -2 s -1 ) 



1.0±0.3 



\EW\ (eV) 



13±5 

i-4±J;g x io 22 

Photoionized Abs. 

Col. dens, (cm* 2 ) 





log 4 (UNITS) 




3.73 ± 0.13 


z (rel. to systemic) 




—0 f|QQ + 0 006 

-u.uyo _ 0 002 


Note. 

details. 


— Spectral parameters for best-fit models which attempt to model the 7.68 keV absorption feature. See text for 


Table 8 


Model 12 fits to each half of the EPIC-pn spectrum 




First half 

Second half 

E 2-10 (erg 

—2 — 1 
cm * s 

9.4xl0 -u 

9.2xl0" u 

X 2 /do f 


1883.15/1887 

1933.52/1887 

r 


1 707+0.031 
-0.013 

1 77 ^+ 0.022 
A -' /D -0.017 

Red peak 

E (keV) 

6-441 ™ 2 

6 45 + 0,02 
D «^°_o.oi 

I ( 10 -5 ph cm -2 s -1 ) 

8-0 ±g;l 

0 q+o.9 
».4_o 7 


EW (eV) 

73l? 

77±8 

Blue peak 

E (keV) 

7 03 +01 ° 
/.uo _ 0 13 

7 03 +0 1 ° 
' • UJ — 0.12 

0 


2 g+0.6 
ZD -1.0 

2 n+0.3 
z - u - 0.6 

Ain 


44±19 

20 ^ 2 j 

Inch 


1Q + 8 O 
AU -10 

q+9o 

U -0 

Gaussian 

E (keV) 

7.68l°°® 

7 cq+0-08 

' -Uo_ 0 .06 


cr 

531 If 

is±ll 2 


| / | ( 10 -5 ph cm -2 s -1 ) 

1 0+O.8 
x ‘°-0.7 



\EW\ (eV) 

16±l° 

10 ±| 


F-test to include 

99.6% 

92.8% 


Note. — Best-fitting model 12, featuring dual-disklines for the Fe K emission peaks and an inverted Gaussian for the 7.68 keV 
absorption feature, to the first and second halves of the EPIC-pn data. All spectral parameters are consistent with remaining 
constant over timescales <136 ksec. The final row shows the level of significance to include the inverted Gaussian according to 
an F-test used in the ’standard’ manner (see text for details). 
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Table 9 

F-TEST RESULTS TO THAW r OR he IN THE TIME-RESOLVED FITS 


Time 

Scale 

r 

F 

r 

Prob 

^Fe 

F 

^Fe 

Prob 

Short 

0.35 

0.95 

0.02 

1.00 

Medium 

5.53 

2.3x 10 -12 

1.51 

7.8xio- 2 

Long 

3.04 

2.3x 10~ 5 

2.28 

1.9xl0~ 3 


Note. — Results of F-tests to determine which parameters it was significant to thaw in the time-resolved fits. High values of 
the F-statistic and low-values of the probability (of observing that value of F from a random set of data) indicate that the fits 
show improvement to thaw that parameter, suggesting statistically significant variations. 


Table 10 

Mean spectral fit variability parameters and errors 


Time 

Scale 

Mean F 2-10 
10 -11 ph cm -2 s -1 

Mean 

r 

Mean he 

10 -5 ph cm -2 s _1 

r 

r 

r 

Prob 

he 

r 

he 

Prob 

Short 

9.91±0.12 

1.822±0.029 

8.30±0.76 

0.803 

1.6xl0~ 2 

-0.106 

0.80 

Medium 

13.29±0.06 

1.828±0.033 

20.47±3.43 

0.819 

3.2x 10 -5 

0.272 

0.27 

Long 

15.62±0.09 

1.917±0.067 

16.60±4.12 

0.656 

3.1xl0- 3 

-0.059 

0.82 


Note. — Mean values of F 2 - 10 , T and he- See text for details on error calculations. The Pearson correlation coefficients r 
in Columns (4) and (6) are for F and he, respectively, plotted against F 2 - 10 , as shown in Figure 11. The corresponding null 
hypothesis probabilities (of achieving that value of r from a random set of data) are listed in Columns (5) and (7). See §4.4 for 
notes relating PCA and pn flux normalizations. 


Table 11 

Fear VALUES 


Time 

Continuum 

Line 

Scale 

Fvar (%) 

Fvar (%) 

Short 

3.50 ± 0.45 

Undef. 

Medium 

9.34 ± 0.14 

3.49 ± 13.99 

Long 

12.99 ± 0.17 

5.54 ± 19.42 


Note. — Undefined fractional variability measurements indicate a measured variance that is smaller than that expected solely 
from measurement noise. 
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Fe K Emission and Absorption in IC 4329a 



Fig. 1. — 10-12 keV pn background light curve. The dotted line at 0.3 count s _1 corresponds to a count rate of B + 2cr, where B is the 
mean background rate and a is the standard deviation of the light curve. We tested for background contamination by filtering out data taken 
when the background count rate exceeded this threshold. However, we found no impact on the Fe K profile, the focus of this paper, and used 
the unfiltered data for maximum signal-to-noise in the analysis. 
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Fig. 2. Panel a) shows the RXTE PCA, HEXTE A cluster and HEXTE B cluster spectra, in black, red and blue, respectively. Panel 
b) shows the data/ model residuals when the data are fit with Model 1, an absorbed power-law model. Panel c) shows the residuals when a 
Gaussian is included to model Fe K emission (Model 2). Panel d) shows the residuals when a PEXRAV reflection component is also included. 
For clarity in panels b), c), and d), residuals in x-space are plotted (typical uncertainties on the HEXTE data/model ratio points were ±0.2 
and so data/model ratios are not plotted). 




ata/model normalized counts 
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Fig. 4. — Data/model residuals in the Fe K bandpass when various models are fit to the EPIC pn data. All plots here denote models 
which include a power-law and a PEXRAV reflection component with R fixed at -0.51, with all components modified by neutral absorption. 
Data have been rebinned by a factor of 15. Panel a) shows the residuals using just the power-law and PEXRAV component (e.g., the ~7.1 keV 
edge has been fit compared to the residuals in Figure 3b). Panel b) shows the residuals to Model 3, which included a Gaussian to model Fe 
Ka emission. Panels c) and d) show the residuals to Models 4 and 5, double Gaussian models for the Fe Ka and K/3 emission lines. The 
energy centroid for the K/3 line was free in Model 4, but fixed at 7.06 keV in Model 5. Panel e) shows the residuals to Model 7, an attempt 
to model both peaks with a single Laor diskline component; there are large residuals suggesting a poor fit to the data. Panel f) shows the 
residuals to Model 8, wherein each peak was modeled with a separate Laor component. Note that the residuals near 5. 8-6.1 keV and 6. 7-6. 9 
keV are improved compared to Model 5, suggesting that the dual-LAOR model is.fitting red wing emission unmodeled by the double-Gaussian 
model (see text for quantification of this statement). Note also that the absorption-like residuals near 7.68 keV are unchanged in panels a)-f), 
suggesting that the absorption feature persists regardless of how the emission profile is modeled. Panel g) shows the residuals to Model 12, 
which features dual-LAOR disklines to model the emission lines and a narrow inverted Gaussian to model the 7.68 keV absorption feature. 
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Fe K Emission and Absorption in IC 4329a 



6 7 8 


Energy (keV) 

Fig. 5. — Derived confidence contours of line intensity versus rest-frame energy when a narrow Gaussian is added to Model 8, the dual-LAOR 
diskline model. The upper panel is for a Gaussian in emission, and suggests there is no additional emission (at 90% confidence or greater) 
to be modeled. The lower panel is for a Gaussian in absorption, and suggest (at >99% confidence) a narrow absorption feature at 7.68 keV. 
Solid, dashed, dashed-dotted, and dotted lines denote 68%, 90%, 95% and 99% confidence levels, respectively. The best-fitting parameters to 
Model 11 are marked with a cross.See text for details. 



Fig. 6. — The pn source (top) and background (bottom) spectra, rebinned by a factor of 20, showing that the source is well over an order 
of magnitude higher in counts s -1 keV -1 compared to the background in the Fe K bandpass. 
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Fig. 8. — Data/model residuals to a simple absorbed power-law (rebinned by a factor of 10) for the MOS2 (black) and pn (grey), showing 
that the MOS2 data are consistent with the pn. 
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Fe K Emission and Absorption in IC 4329a 



Energy (keV) 

Fig. 9. — Data/model residuals for the first half (black) and second half (gray) of the pn data when fitting to a model consisting of a 
power-law and a PEXRAV reflection component, all modified by neutral absorption. All derived spectral parameters are consistent with each 
other, suggesting no evidence that the profile varies on short time scales. 
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Fig. 10. — Light curves for the 2-10 keV continuum flux F 2-10 (top), photon index T (middle), and Fe K flux Ip e (bottom) for all three 
time scales, derived from time-resolved spectroscopy. Fh-io is ' n units of 10~ n erg cm -2 s _1 . / Fe > s in units of 10 -s ph cm -2 s -1 . 
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Fig. 11. — Photon index T (top) or Fe K flux /p e (bottom) plotted against 2-10 keV continuum flux F 2-10 for all three time scales, derived 
from time-resolved spectroscopy. F 2-10 is in units of 10 -u erg cm -2 s -1 . /Fe is in units of 10 -5 ph cm -2 s _1 . These plots, along with 
Figure 10 and Tables 10 and 11, suggest that temporal variations in F are well-correlated with those in F 2 _io on all three time scaled studied, 
but there is not strong evidence for Ip e to vary on any of the time scales studied. See text for details. 





